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1.0  Introduction 


Aerospace  vehicles  of  the  present  and  future  are  being  designed  to 
operate  in  extremely  large  flight  envelopes  in  terms  of  the  velocity, 
acceleration,  altitude,  and  maneuverability  that  they  are  required  to  achieve. 
Not  only  are  these  vehicles  required  to  operate  in  such  large  flight  envelopes, 
but  they  must  maintain  specified  performance  margins  while  doing  so.  The 
development  of  such  vehicles  presents  a  great  challenge  to  the  flight  control 
system  designer.  The  lack  of  adequate  analytical  techniques  and 
experimental  data  for  systems  expected  to  operate  on  the  extremities  of 
existing  operational  boundaries  results  in  significant  uncertainty  concerning 
the  dynamical  description  of  the  vehicle.  Furthermore,  aeroservoelastic 
effects  arising  from  flexible  vehicle  characteristics  will  become  more 
pronounced  and  unpredictable  at  higher  velocities  and  accelerations.  This 
means  that  the  mathematical  models  currently  used  to  describe  vehicle 
motion  and  to  design  flight  control  systems  may  not  be  appropriate  in  terms  of 
parametric  and  structural  fidelity  for  many  of  the  operating  conditions  of 
interest. 

Adaptive  feedback  control  has  long  been  regarded  as  one  of  the  leading 
choices  for  designing  control  laws  for  systems  that  are  subject  to  the  types  of 
uncertainties  described  above.  However,  most  of  the  existing  adaptive 
control  techniques  have  performance  limitations  and/or  implementation 
drawbacks  that  forbid  them  from  providing  the  operating  robustness  and/or 
operational  confidence  required  of  flight  control  systems  for  realistic  vehicles. 
The  primary  goal  of  this  Phase  I  effort  was  to  develop  a  control  system  design 
methodology  that  would  produce  a  control  law  capable  of  providing  the 
necessary  performance  in  the  presence  of  uncertainties  while  avoiding  the 
limitations  and  drawbacks  of  current  adaptive  control  design  techniques. 
Dynetics'  innovative  approach  was  to  apply  Disturbance  Accommodating 
Control  (DAC)  thor>r\y  to  the  flight  control  pi-oblem  to  create  a  novel  alternative 
to  traditional  adaptive  control  system  design. 

Controllers  that  were  designed  using  this  innovative  approach  were 
compared  to  other  controllers  designed  using  various  contemporary  adaptive 
control  theories.  These  designs  were  based  upon  Self-Tuning  Regulator 
(STR)  and  Model  Reference  Adaptive  Control  (MRAC)  concepts.  All 
controller  designs  were  evaluated  in  identical  operating  environments  which 
included  parameter  variations,  unmodeled  dynamics,  and  external 
disturbances.  The  relative  merits  of  each  technique  were  compared  in  terms 
of  performance  robustness,  design  complexity,  and  operational  confidence. 

In  every  case  the  controllers  designed  using  Dynetics’  innovative  approach 
were  able  to  equal  or  surpass  the  STR  and  MRAC  controllers  in  terms  of 
performance  robustness  while  preserving  a  linear  time-invariant 
implementation  structure. 
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The  development  of  a  design  methodology  that  will  produce  realizable 
control  laws  that  can  compensate  for  significant  uncertainties  is  a  contribution 
to  the  field  of  teedoack  control  theory  in  general,  and  to  the  aerospace 
industry  in  particular.  The  feasibility  of  using  DAC  principles  to  develop  such 
an  adaptive  control  system  design  methodology  was  demonstrated  with 
complete  success  in  this  Phase  I  effort.  The  applications  of  the  resulting 
design  methodology  are  as  varied  as  the  types  of  systems  under 
development  and  consideration.  All  practical  control  systems  are  designed 
by  using  a  mathematical  model  that  has  been  simplified  to  a  degree.  Even  in 
the  most  comprehensive  cases  there  are  invariably  higher  order  effects  that 
are  deemed  to  be  negligible  in  this  sense,  all  controller  designs  are 
subjected  to  varying  degrees  of  uncertainty  when  implemented  in  reali.-^tic 
systems  and  therefore  will  benefit  from  the  results  of  this  study. 

Evaluation  of  design  complexity,  confidence  of  operation  and 
performance  robustness  verifies  that  the  DAC  based  designs  are  superior  to 
the  other  adaptive  control  techniques.  The  MRAC  system  exhibited  very 
good  performance  characteristics  for  certain  choices  of  tuning  parameters, 
but  the  inherent  nonlinearity  of  the  resulting  controller  tends  to  limit  the 
applicability  of  this  method.  The  only  tools  available  to  analyze  and  measure 
the  reliability  and  robustness  of  the  closed  loop  system  are  repeated 
simulation  and  Liapanov  stability  analysis.  Again,  the  STR  approach 
resulted  in  a  nonlinear  closed  loop  system  which  was  subject  to  the  same 
analysis  problems  as  the  MRAC.  Beyond  the  issue  of  nonlinearity  is  the 
fundamental  limitation  of  the  estimation  of  system  parameters  in  the  presence 
of  external  disturbances.  The  resulting  DAC  based  controllers  are  linear 
time-invariant  systems,  which  means  that  all  of  the  standard  linear  analysis 
tools  (Nyquist,  Bode,  root  locus,  etc.)  are  still  applicable. 


2.0  Objectives 


The  objective'^  of  the  Phase  I  study  were  as  follows: 

(1 )  Chafr.wtenzation  of  parameter  and  dynamical  uncertainties 

(2)  Application  of  Disturbance  Accommodating  Control  theory  to 
dovelopmeni  of  control  law  design  methodologies  for  systems  subject  to 
Significant  uncertainties 

(3)  Specification  of  candidate  control  algorithms 

(4)  Problem  formulation  and  candidate  solution  development,  and 

(5)  Performance  demonstration. 

The  nature  of  the  uncertainties  that  a  control  system  will  encounter  can 
be  characterized  by  considering  three  types  of  effects;  parameter  uncertainty, 
uncertain  external  disturbances,  and  unmodeled  dynamics, 

2.1  Characterization  of  Uncertainties 

Parameter  uncertainty  can  be  mathematically  described  by 
consideration  of  the  general  state-space  system  equation  shown  in  Figure 
2.1  Parameter  uncedamty  is  manifested  by  perturbations  away  from  some 
expected  (nominal)  values  of  the  elements  of  the  A,  B,  and  C  system  matrices 
These  perturbations  AA,  AB,  and  AC  also  correspond  to  the  uncertainty 
associated  with  the  coefficients  of  the  differential  equations  which  govern  the 
behavior  of  the  system  when  expressed  in  a  transfer  function  form  as  shown 
in  Figure  2.2. 

it  IS  interesting  to  note  the  correspondence  between  the  state-space 
system  representation  and  the  transfer  function  representation  in  terms  of  the 
numerator  and  denominator  coefficients  and  their  respective  locations  in  the 
A.B.C  system  matrices  The  transfer  function  of  Figure  2.2  can  be  written  m 
several  different  canonical  state-space  forms,  each  of  which  has  a  distinct 
configuration.  In  most  cases  the  denominator  coefficients  will  appear  as 
elements  of  the  A  matrix  (either  in  a  single  row,  a  single  column,  or  some 
combination  of  each).  The  numerator  coefficients,  however,  can  appear  as 
entries  in  the  C  matrix  alone,  in  the  B  matrix  alone,  or  as  some  combination  of 
the  two. 

On  the  surface  it  appears  that  the  choice  of  a  canonical  form  is 
completely  arbitrary  and  theretore  should  not  affect  the  controller  design 
process.  This  is  in  fact  the  case  when  all  of  the  parameters  of  the  system  are 
certain  or  known.  However,  if  the  parameters  are  subject  to  perturbations  or 
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i  =  Ax  +  Bu 


y  =  Cx 


Figure  2.1.  Parameter  Uncertainty,  State-Space  Representation 


Numerator  Uncertainty 


Denominator  Uncertainty 


Figure  2.2.  Parameter  Uncertainty,  Transfer  Function 
Representation 
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uncertainties,  it  is  very  important  to  understand  how  the  uncertainty  will  affect 
the  system  response.  This  in  turn  dictates  that  the  actual  state-space 
description  of  the  system  should  have  some  relative  significance  and 
therefore  should  not  be  chosen  arbitrarily.  More  specifically,  the  B-matrix 
describes  a  mapping  of  the  system  input  into  the  particular  state-space 
described  by  the  A-matrix,  in  a  geometrical  sense.  And  conversely,  the  C- 
matrix  describes  the  mapping  of  the  entire  state  vector  into  some  observed 
sub-space.  Or  stated  more  heuristically,  the  B-matrix  describes  how  the  input 
gets  distributed  into  the  system,  and  the  C-matrix  describes  how  the  system 
motion  gets  described  in  the  form  of  observed  outputs. 

Another  type  of  uncertainty  is  external  disturbances.  External 
disturbances  are  defined  in  the  context  of  this  study  as  those  inputs  which  are 
both  directly  unmeasurable  and  uncontrollable  by  the  control  system 
designer.  External  disturbances  are  mathematically  described  by  the 
addition  of  forcing  terms  in  the  system  equations  as  shown  in  Figure  2  3. 
These  disturbances  can  also  be  described  in  transfer  function  or  blocK 
diagram  form  by  the  configuration  shown  in  Figure  2.4. 

The  third  type  of  uncertainty  considered  in  this  study  was  unmodeled 
dynamics.  Unmodeled  dynamics  can  be  described  mathematically  by  the 
system  equations  shown  in  Figure  2.5.  The  resulting  effect  of  the  unmodeled 
dynamics  is  that  the  dimension  of  the  system  model  used  to  design  the 
control  system  is  less  than  the  dimension  of  the  actual  system  to  be 
controlled.  Unmodeled  dynamics  can  be  further  divided  into  two  categories: 
additive  and  multiplicative,  as  shown  in  Figure  2.6.  The  effect  in  both  cases  is 
similar  in  that  each  type  increases  the  dimension  of  the  actual  system,  but 
their  treatment  in  the  design  of  a  control  system  is  very  different. 

2.2  General  Problem  Statement 

One  of  the  original  motivations  for  this  study  was  to  develop  a  ro-^Tol  law 
design  methodology  for  hypervelocity  vehicles,  therefore  the  class  of 
problems  chosen  for  consideration  was  taken  to  be  representative  of  these 
types  of  aerospace  vehicles.  A  benchmark  design  problem  was  formulated  to 
capture  the  important  features  of  designing  a  realistic  flight  control  system  for 
such  a  vehicle. 

The  particular  problem  formulated  was  to  design  a  controller  for  the 
longi\  :  .lal  motion  of  an  aerospace  vehicle.  This  problem  was  used  as  the 
bas:  T  test  problem  throughout  this  Phase  I  study.  A  linearized  dynamical 
mocei  aescribing  the  vehicle's  motion  in  terms  of  angular  rate  about  some 
center  ;  avity  reference  was  developed.  The  transfer  function  which 
ch.  .rrctCiizes  this  dynamical  model  is  shown  in  Figure  2.7.  This  dynamical 
model  is  very  common  and  is  found  to  reoccur  frequently  m  a  multitude  of 
aerospace  control  problems.  This  particular  model  describes  the  pitch  or  yaw 
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X  =  Ax  +  Bu  +  Fw  +  . . . 
y  =  Cx  +  Hw  +  ... 
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Figure  2.3.  External  Disturbances,  State*Space  Format 


disturbance 


Figure  2.4.  External  Disturbances,  Transfer  Function  Format 
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Figure  2.5.  Unmodeled  Dynamics 


•  VEHICLE  PITCH  PLANE  DYNAMICS 
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Figure  2.7.  Dynamical  Description  of  Baseline  Problem 
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rate  for  aircraft,  missile  angular  turning  rate  (either  tail  or  canard  controlled), 
and  space  launch  vehicle  angular  turning  rate  (atmospheric  portion  of  flight). 
The  use  of  this  particular  dynamical  problem  makes  the  results  of  this  Phase  I 
study  very  general  and  applicable  to  a  variety  of  systems. 

The  benchmark  problem  was  to  design  a  suitable  controller  which  could 
maintain  a  set  of  specified  performance  requirements  in  the  presence  of  the 
uncertainties  described  in  section  2.1. 
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3.0  Performance  versus  Stability 

A  succinct  distinction  is  made  between  performance  and  stability  for  the 
controller  designs  considered  in  this  study.  For  our  purposes  stability  is 
defined  in  the  traditional  manner  by  describing  the  controllers  ability  to  drive 
the  homogeneous  error  system  to  zero  asymptotically.  Performance  is 
defined  in  terms  of  the  controllers  ability  to  force  the  system  output  to 
■  eproduce  a  desired  transient  shape  and  to  maintain  specified  input 
command  tracking  margins. 

This  distinction  between  stability  and  performance  is  necessary  due  to 
the  fact  that  almost  all  of  the  research  in  adaptive  control  for  the  last  three 
decades  has  been  preoccupied  with  proving  global  stability  for  certain 
classes  of  plant  uncertainty.  Engineers  concerned  with  the  application  of 
adaptive  control  techniques  are  much  more  concerned  with  the  controllers 
ability  to  maintain  performance  specifications  in  the  presence  of  uncertainties 
which  may  be  unknown  exactly  but  which  are  invariably  bounded  in  realistic 
problems.  For  this  reason,  the  benefits  and  limitations  of  various  adaptive 
control  designs  may  not  be  at  all  apparent  to  the  flight  control  system 
designer. 

In  the  area  of  flight  control,  performance  is  typically  the  ultimate 
justification  for  using  any  type  of  active  feedback  control.  The  airframes 
encountered  are  almost  invariably  statically  stable,  and  only  in  very  recent 
history  have  there  been  any  actual  systems  built  which  rely  on  active  stability 
augmentation  for  operation.  Active  control  is  employed  on  virtually  all 
vehicles  as  a  means  of  enhancing  the  performance  of  the  closed  loop  system. 
In  the  case  of  piloted  aircraft,  performance  translates  into  the  specific 
handling  qualities  of  the  vehicle  which  are  required  by  the  pilot.  For  a  missile, 
the  flight  control  system  performance  can  be  directly  measured  in  terms  of  the 
control  system  contribution  to  overall  miss  distance.  For  a  space  launch 
vehicle,  performance  is  specified  in  terms  of  the  control  systems  ability  to 
maintain  robustness  margins  due  to  the  fundamental  importance  of  achieving 
mission  success.  In  each  of  these  realistic  cases,  it  is  obvious  that  stability  is 
assumed  and  that  performance  is  really  the  primary  concern. 

An  example  of  a  specific  performance  template  which  might  describe  the 
required  closed-loop  pitch-rate  response  for  a  piloted  aircraft  is  shown  in 
Figure  3.1 .  Both  the  transient  and  tracking  requirements  are  apparent  from 
the  specified  regions  of  acceptable  and  unacceptable  performance.  This 
specification  is  typical  of  the  type  of  handling  quality  requirement  available  in 
several  Military  Standard  documents  for  piloted  aircraft. 
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Figure  3.1.  Typical  Performance  Boundaries 
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4.0  Disturbance  Accommodating  Control 

Disturbance  Accommodating  Control  (DAC)  theory  has  existed  as  a 
design  discipline  since  the  late  1960's  [1].  Concepts  and  techniques  for  the 
characterization  and  design  of  disturbance  accommodating  controllers  have 
been  refined  and  extended  into  a  fairly  complete  technology  [2],  [3],  [4], 

The  original  motivation  for  the  development  of  DAC  was  to  generalize 
the  disturbance  compensation  techniques  of  classical  control,  such  as 
integral  feedback,  feedforward  compensation,  and  notch  filtering,  into  the 
framework  of  multivariable  systems.  Figure  4.1  illustrates  a  typical  closed 
loop  system  consisting  of  a  controller  whose  purpose  is  to  force  the  plant  to 
follow  a  servo-command  in  the  presence  of  plant  uncertainty,  external 
disturbances,  and  corrupt  measurements.  Disturbances  are  defined  as  all  of 
those  plant  and  controller  external  inputs  or  internal  actions  that  are  totally 
unexpected  and  uncontrollable  from  the  designers  point  of  view. 

Disturbances  can  be  generally  classified  into  two  broad  categories: 
erratic  (noise-like)  disturbances  and  waveform-like  disturbances  (See  Figure 

4.2. )  Erratic  noise  signals  are  best  described  using  statistical  properties  such 
as  mean  and  variance.  Waveform  type  disturbances  are  typically  non- 
ergodic  in  nature  and  therefore  their  statistical  properties  are  not  useful  in 
terms  of  designing  a  compensating  control  strategy.  A  specific  example  of 
this  is  the  case  of  a  constant  random  disturbance.  If  an  ensemble  average  is 
constructed,  the  mean  value  could  well  be  zero,  but  if  a  temporal  average  is 
calculated  the  mean  value  is  simply  the  value  of  the  disturbance  itself.  By 
taking  a  temporal  average  the  mean  value  with  respect  to  time  is  calculated, 
which  in  the  case  of  a  random  constant  is  simply  the  value  of  the  constant.  An 
ensemble  average  is  calculated  by  taking  a  cross-section  of  the  entire 
possible  sei  of  random  constants  at  a  particular  instance  in  time,  which  gives 
the  mean  value  of  the  sample  set  as  opposed  to  the  mean  value  of  a 
particular  sample.  In  this  case  knowledge  of  the  statistical  properties  of  the 
disturbance  are  not  useful  for  designing  a  control  strategy  to  offset  the  effects 
of  the  constant  disturbance.  In  fact  it  is  well  known  that  integral  feedback  is 
the  proper  compensation  scheme  for  unknown  constant  random 
disturbances. 

The  central  theme  of  DAC  theory  is  in  modeling  disturbance  processes 
as  time  signals  which  have  particular  "waveform"  characteristics  (See  Figure 

4.3. )  These  types  of  signals  can  always  be  mathematically  described  by  the 
following  type  of  equation 


co{t)=c,f,(t)  +  C2f2(t)  + . +  c^U^) 


This  is  a  function-space  representation  of 
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PLANT 

UNCERTAINTY 


Figure  4.1.  Block  Diagram  of  Closed  Loop  System 
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to(t)  where  the  basis  set  {f, ,  fp-  •  ■  ^m)  'S  chosen  to  reflect  the  natural 
waveform  patterns  inherent  in  the  disturbance  behavior.  These  functions  fj 

embody  the  character  of  the  disturbance  signal.  The  choice  of  the  particular 
types  of  functions  to  be  used  as  the  basis  set  is  strictly  a  design  alternative. 
These  functions  could  be  constants,  sinusoids,  exponentials,  polynomials, 
etc.  The  Cj's  are  constants  which  are  chosen  once  the  basis  set  is  defined. 

The  particular  values  used  for  these  constants  specifies  the  exact  time 
characteristic  of  the  disturbance  model.  This  is,  in  effect,  a  generalized 
spline-iunction  model  of  the  disturbance. 

Although  DAC  was  initially  developed  as  a  means  of  compensating  for 
external  disturbances,  it  has  also  been  suggested  as  a  means  of 
compensating  for  other  system  uncertainties  such  as  parameter 
perturbations.  Consider  the  following  dynamical  system 


X  =  Ax  +  Bu  +  Fw 


y  =  Cx 

Parameter  uncertainty  is  represented  by  perturbations  in  the  A,  B,  and  C 
system  matrices,  away  from  their  nominal  or  expected  values.  If  the  nominal 

values  are  written  as  A^,,  B^,,  and  and  the  perturbations  are  5A,  5B,  and  6C. 
then  the  system  equations  can  be  rewritten  as 

X  =  (An  +  5A)x  +  (Bn  +  8B)u  +  F.v 
y  =  (Cn  5C)x 


and  the  uncertainties  can  be  grouped  as 

X  =  AnX  +  BnU  +  (8Ax)  +  (8Bu)  +  Fw 
y  =  CnX  +  (8Cx) 


The  basic  philosophy  of  the  DAC  compensation  technique  is  to  treat  all  of  the 
uncertainties  shown  in  the  system  equations  above  as  some  type  of 
undesirable  disturbance  (either  internal  or  external).  Following  standard 
DAC  design  procedures,  these  uncertainties  are  treated  as  the  output  of  a 
dynamical  process.  The  next  step  in  the  design  procedure  is  to  construct  a 
state-space  model  of  the  disturbance  process. 
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Let 


W  —  +  C2f2  +  ••■  +  Cp-iffYi 


In  standard  DAC  practice  the  disturbance  model  is  constructed  by  first 
choosing  a  set  of  basis  functions  for  the  spline  model  which  characterize  w. 

In  most  cases  the  designer  will  have  little  if  any  knowledge  concerning  the 
waveform  characteristics  of  the  process  w,  and  therefore  would  have  no 
rational  way  to  assign  a  set  of  effective  basis  functions  to  the  spline  model.  In 
this  case,  practice  has  shown  it  to  be  very  effective  to  use  polynomials  (linear, 
quadratic,  cubic,  etc.)  for  the  basis  functions.  This  follows  from  similar 
practices  in  the  area  of  approximation  theory.  Quadratic  and  cubic  splines 
can  be  fit  to  waveform  types  of  signals  to  high  degrees  of  accuracy  over 
specified  intervals  of  interest. 

The  state-space  model  for  w  is  constructed  by  solving  the  inverse 
problem  from  differential  equation  theory,  i.e.  find  the  differential  equation  to 
which  w  is  a  solution.  Once  the  differential  equation  in  w  is  formulated,  then 
the  state-space  model  of  this  differential  equation  represents  the  state  model 
of  the  disturbance  process.  As  an  example  consider  a  case  when  very  little  is 
known  about  the  actual  uncertainty  and  the  designer  chooses  to  employ  a 
quadratic  spline  to  characterize  w.  Then 

w=c,  +C2t  +C3t^ 

which  IS  the  general  solution  to  the  following  differential  equation 

w  =  0 

which  can  be  represented  by  the  following  state-space  model 
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w  =  [l  0  O] 


Or,  rewritten  in  general  matrix  format 

z  =  Dz 
w  =  Hz 

In  the  case  that  nothing  at  all  is  known  about  the  uncertainty,  the  polynomial 
spline  approach  is  always  a  valid  design  solution.  But  for  the  case  of  5A 
perturbations  in  particular,  we  see  that  the  natural  choice  for  the  basis 
functions  of  the  spline  are  simply  the  eigenfunctions  of  the  nominal  system  [4). 
This  choice  of  basis  functions  has  been  called  the  Adaptive  DAC  (ADAC) 
approach.  The  distinction  between  choices  of  basis  functions  and  the  design 
alternatives  available  are  shown  graphically  in  Figure  4.4. 

In  either  case  the  design  proceeds  in  similar  fashion  once  the  basis 
functions  have  been  chosen.  The  basic  idea  in  compensating  for  the 
uncertainty  characterized  by  z  is  to  identify  it’s  effect  and  to  cancel  it  (there  are 
other  possible  strategies  such  as  minimization  of  uncertainty  or  possible 
utilization  of  uncertainty  [2],  but  only  cancellation  will  be  considered  here).  To 
cancel  the  effects  of  z  it  must  first  be  identified  in  some  manner  In  DAC  this 
identification  is  achieved  by  means  of  a  composite  state  observer  constructed 
by  augmenting  the  original  system  with  the  dynamical  model  of  z 
Specifically,  recall 


X  =  A^x  +  Bu  +  Fw 
y  =  Cx 


and 


z  =  Dz 


w  =  Hz 


then 
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-  DERIVE  "WAVEFORM  MODEL"  OF  o){t) 

oo(t)  =  (t)  +  C2  f2(t)  +  ■■•  + 
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(JUMP  ONCE-IN-A-WHILE) 

KNOWN  (CHOSEN)  SET 
OF  "BASIS-FUNCTIONS" 


Figure  4.3.  Disturbance  With  Waveform  Characteristics 


disturbance  accommodating  control 

a)(t)  =  f^(t)  *  C2  ••• 


CONVENTIONAL  DAC 

•  BASIS  FUNCTIONS  OF  SPLINE  ARE 
CHOSEN  BASED  ON  EXPECTED 
CHARACTERISTICS  OF  DISTURBANCE 

•  IF  THE  NATURE  OF  THE  DISTURBANCE 
IS  UNKNOWN.  POLYNOMIALS  (LINEAR, 
QUADRATIC,  CUBIC)  PERFORM  NICELY 


ADAPTIVE  DAC 

•  BASIS  FUNCTIONS  OF  SPLINE 
ARE  CHOSEN  AS  THE  -NATURAL- 
EIGENFUNCTIONS  OF  THE 
NOMINAL  SYSTEM 

•  THIS  IS  THEORETICALLY  DERIVED 
FOR  THE  CASE  OF  PARAMETER 
PERTURBATIONS  ONLY 


Figure  4.4.  Disturbance  Accommodating  Control  Design 

Alternatives 
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This  state-space  system  is  simply  another  linear  system  for  which  a 
composite  observer  can  be  generated.  The  composite  observer  can  be 
written  as 
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A  block  diagram  of  the  composite  observer  is  shown  in  Figure  4.5. 

4.1  Adaptive  DAC  Design  for  Baseline  Problem 

For  specific  examples  of  the  design  procedure,  again  consider  the 
benchmark  problem  outlined  in  section  2.2.  The  plant  is  given  by 

y _ ks  +  b 

u  “  2 

s  +  825  +  a^ 


which  corresponds  to  the  following  differential  equation 


y  +  +  a,y  =  ku  +  bu 


This  system  can  be  written  in  the  following  state-space  format 


'  0  1  ' 

k 

- 1 

CM 

X 

1 

-a  -a 

L  1  aJ 

I - 

X 

+ 

a^k  -  b 

L  2  J 

Now  assume  that  the  objective  is  to  track  some  constant  input  y^p.  The 
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Figure  4.5.  Block  Diagram  of  Composite  State  Observer 


following  error  can  then  be  defined 


f-ysp-y 

y=ysp-f 

and  the  following  error  system  is  written 
-e-  a^-  a^e+a^ysp  =  ku +bu 

e  +  a^e  +  a^e  =  a^Ysp  -  ku  -  bu 

If  the  plant  parameters  are  perturbed  from  their  nominal  values 

a,  =(a,„  +  Sa,) 

then  the  error  system  becomes 

£  +  a^  +  a,e  =  a,ysp-5a2e-  2a,e-ku-bu 


Now  let 


z^  =  -(ba^z+  6a, e) 

=  ku  +  (any  other  disturbance  terms) 

where  Zg  represents  an  effective  “disturbance"  due  only  to  the  effects  of 

parameter  variations,  and  z^j  represents  a  disturbance  which  contains  all 

other  external  and  internal  uncertainty  effects, 
then 


e+a2e>a,e  =  a,y3p-bu  +  23-z^ 

Recall  that  the  objective  is  to  force  y  to  follow  a  constant  input,  or  to  force  e  to 
zero  in  a  prescribed  manner.  Examination  of  the  error  system  equation  as 
written  above  indicates  exactly  the  control  law  necessary  to  force  e  to  zero, 
namely 
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The  control  law  is  separated  into  three  parts  Up,  Ug,  and  u^j.  The  primary 
control  signal  Up  is  designed  to  achieve  the  fundamental  objective  of  tracking 

the  set-point  command.  This  part  of  the  control  signal  is  designed  assuming 
that  the  plant  is  free  of  disturbances,  and  is  chosen  exactly  the  same  as  in 
other  state-space  based  design  approaches  such  as  pole  placement.  In  this 
particular  problem  no  pole  placement  was  necessary  due  to  the  fact  that  the 
nominal  system  characteristic  equation  is  exactly  the  same  as  the  desired 
closed-loop  characteristic  equation.  The  signal  Ug  is  chosen  to  cancel  the 

effects  of  parameter  variations  in  the  closed-loop  system.  The  signal  up  is 

designed  to  cancel  the  effects  of  all  other  internal  and  external  uncertainties. 
Due  to  the  fact  that  the  composite  observer  and  therefore  the  resulting 
compensator  are  linear  systems,  the  property  of  superposition  is  valid  which 
allows  the  design  of  the  control  law  to  be  separated  into  three  parts  which  can 
be  performed  independently.  One  of  the  fundamental  attributes  of  DAC 
design  is  the  explicit  identification  of  quantities  such  as  zp  and  Zg  which 

facilitates  the  type  of  control  law  design  shown  above. 

Now  all  that  is  needed  to  implement  this  control  law  are  the  estimates  of 
Z3  and  z^.  For  this  the  composite  state  observer  is  necessary.  First  the  state 
models  of  the  uncertainties  must  be  constructed.  Since  the  term  zp 

represents  all  of  the  uncertainty  effects  other  than  parameter  variations,  the 
specific  character  of  this  signal  is  unknown.  Therefore,  we  will  follow 
conventional  DAC  techniques  discussed  earlier  and  will  choose  a  linear 
polynomial  spline  model  for  z^.  This  choice  of  basis  functions  is  a  very  good 
general  model  for  completely  unknown  types  of  signals.  The  corresponding 
state  model  is 
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Now  following  the  ADAC  procedure,  assume  that  the  perturbation  effects  are 
character.zed  by 


Za  +  3  +  3  .Zjv  = 
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The  choice  of  the  nominal  system  parameters  as  the  appropriate  estimation 
matrix  in  the  equations  above  follows  from  the  development  in  [4],  where  it 
was  suggested  that  the  evolution  of  parameter  variation  effects  on  the 
nominal  system  could  best  be  characterized  by  the  eigenfunctions  of  the 
original  nominal  system.  Now  the  complete  composite  system  is  formed 
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An  observer  is  constructed  for  this  composite  system  which  is  given  again  by 


X  =  +  Bu  +  kjjfy  -  Cx] 


where 
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This  observer  will  produce  the  estimates  of  Zg  and  necessary  to  implement 
the  control  law,  which  in  its  final  form  is 


u  =  Up  +  ug  +  u,  =-^ysp  +  (^)  Zai  -  ^d1 


4.2  Conventional  DAC  Design 

A  conventional  DAC  controller  design  can  be  constructed  using  a  similar 
procedure  to  that  outlined  above.  As  discussed  earlier,  the  distinction 
between  the  conventional  DAC  and  Adaptive  DAC  designs  lies  with  the 
choice  of  basis  functions  used  for  the  spline  model  of  the  uncertainty.  In  order 
to  better  understand  the  complete  capabilities  of  the  ADAC  approach,  a 
comparison  with  conventional  DAC  design  methods  was  deemed  to  be 
useful. 
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Again  consider  the  plant  described  in  the  baseline  problem 


^  _  ks  -t-  b 
u  ~  2 

S  +  325  +  a, 


which  can  be  represented  by  the  following  state-space  system 


Now  assume  that  all  of  the  uncertainties  can  be  lumped  into  one  term  z, 
which  can  be  characterized  by  the  following  cubic  polynomial 

z  =  c,  +C2t+C3t^  +  c/ 

In  this  case  we  have  assumed  that  all  uncertainty  effects,  including  those  due 
to  parameter  variations,  can  be  represented  by  the  term  z.  Since  this  term 
includes  effects  from  several  different  internal  and  external  sources,  it  is 
assumed  that  little  is  known  about  the  specific  character  of  this  signal. 
Therefore,  as  discussed  earlier,  a  good  choice  for  the  basis  functions  of  this 
disturbance  model  is  a  polynomial.  This  polynomial  satisfies  the  following 
general  differential  equation 

z  =  0 

which  has  the  state-space  model 
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Recall  the  composite  observer  for  the  system  is 

-k=  ylx  +  Bu  +  kJy-C^l 


where 
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c  =  [^  k,  0  0  0  Oj 

The  control  law  for  this  case  is  simply 

u  =  (^)  ysp  ■  z, 


26 


5.0  Self-Tuning  Regulators 

The  concept  of  self-tuning  regulator  adaptive  controller  design  is  much 
as  the  title  would  indicate.  First  an  appropriate  dynamical  model  for  the 
actual  system  is  developed,  secondly  some  control  law  design  algorithm  such 
as  pole  placement  or  linear-quadratic-regulation  (LQR)  is  chosen,  and  lastly 
the  parameters  of  the  dynamical  model  are  estimated  and  used  to  solve  for 
the  required  closed  loop  control  law. 

5.1  Regulator  Design 

There  are  several  alternatives  available  for  designing  the  regulator 
portion  of  the  STR  controller.  One  particular  technique  which  is  very 
straightforward  is  the  polynomial  approach  to  pole  placement.  Given  the 
plant 


A(s)y  =  B(s)u 

where  A(s),  B(s)  are  polynomials  in  s,  y  is  the  plant  output,  and  u  is  the  plant 
input.  The  general  control  equation  can  be  written  as 

L(s)u  =  -P(s)y  +  H(s)yref  ■ 

Substitution  of  the  first  equation  into  the  second  equation  above  results  in 

[A(s)L(s)  +  B(s)P(s))y  =  B(s)H(s)yref 
which  can  be  written  in  transfer  function  form  as 


_y _ B(s)H(s) 

y,e,  ~  A{s)L(s)  +  B(s)P(s) 
For  the  benchmark  problem  we  have 

A(s)  =  +  a^s  +  a, 

B(s)  =ks  +  b 

and  for  the  regulator  design  we  choose 

L(s)  =  s  + 

P(s)  =  k^s  +  k3 

Then 
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(A(s)L(s)  +  B(s)P(s)]  =  (s^  +  823  +  a^)(s  +  k^)  +  (ks  +  b)(k2S  +  k^) 

=  S  +  (82  k2k)S^  +  (8^  +  ^1^2  ^  ^3^  *  k2b)S  +  ■*■  ^3) 


The  next  step  in  the  design  process  is  to  specify  a  characteristic  polynomial 
which  exhibits  the  transient  behavior  desired  from  the  cLsed  loop  system. 
This  characteristic  polynomial  has  the  following  form 


Now  we  equate  the  coefficients  of  the  last  two  equations  above  and  solve  for 
the  controller  gains.  The  result  is  written  as 


^  =  1(2^2  '  ^  ^(^2d  ■ 

k  =  [a  b(a  -  a  )  +  a  k{a,  -  a  )  +  a  (b  -  ka  )]/A 

3^12  3a  1  2d  1  Id  2' 

A  =b^  -  k(a2b  -  ka^) 


The  final  step  is  to  estimate  the  plant  parameters  a-] ,  a2,  k,  and  b,  choose  the 
values  ^2^,  a-]  ^j,  solve  for  k-j ,  k2,  k3  and  the  problem  is  solved. 

This  specific  procedure  outlines  the  polynomial  approach  to  pole 
placement,  but  this  is  by  no  means  the  only  option  available  to  the  designer. 
Other  popular  design  algorithms  such  as  LQR  could  just  as  easily  be  used  in 
this  kind  of  problem. 

5.2  Estimation 

To  implement  the  controller  of  section  5.1.  we  must  obtain  the  estimated 
values  of  the  plant.  Most  of  the  current  literature  on  self-tuning  regulator 
design  encourages  the  separation  and  optimization  of  the  regulator  and 
estimator  functions  individually.  Several  alternative  algorithms  are  available 
for  estimation  which  are  found  in  references  such  as  [5]  and  [6], 

Most  estimation  algorithms  are  based  upon  minimizing  some  function 
which  measures  the  suitability  of  the  estimate  based  on  observed  data.  The 
two  common  measures  of  suitability  are  the  square  of  the  error  between  the 
estimate  and  the  actual  value  which  results  in  the  class  of  least  squares  (LS) 
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estimators,  and  the  maximum  likelihood  (ML)  class  of  estimators  generated 
by  the  minimization  of  a  probabilistic  likelihood  function.  For  all  practical 
considerations  the  difference  between  the  resulting  algorithms  for  both  types 
of  measure  functions  is  negligible. 

Estimation  can  be  performed  in  two  basic  modes,  batch  and  recursive. 

In  the  batch  mode  a  certain  amount  of  input  and  output  data  is  collected  and 
processed  in  an  "off-line"  manner  to  generate  the  required  estimates.  In  this 
mode  it  is  possible  to  perform  the  estimation  in  either  the  time  domain  or  in 
the  frequency  domain.  The  second  mode  of  estimation  is  the  recursive  mode. 
When  data  is  processed  recursively  a  new  estimate  is  produced  each  time  a 
successive  sample  of  input  and  output  data  is  introduced  into  the  estimator. 
The  recursive  mode  can  be  thought  of  as  an  "on-line"  mode  of  operation.  For 
the  type  of  systems  considered  in  this  study  it  is  necessary  that  the  estimates 
be  available  in  an  "on-line"  manner  since  the  nature  of  flight  control  problems 
IS  that  the  system  under  control  is  constantly  changing.  Therefore  it  is 
computationaly  impossible  to  perform  any  kind  of  batch  estimation  using 
typical  avionics  hardware  systems.  For  this  reason  only  recursive  estimation 
algorithms  were  considered  in  this  study. 

Several  commonly  available  Recursive  Least  Squares  (RLS)  estimation 
algorithms  were  evaluated  [7],  [8],  [9].  The  nature  of  these  algorithms  makes 
them  particularly  amenable  to  implementation  as  discrete-time  difference 
equations.  This  necessitates  the  use  of  a  discrete-time  model  for  the  plant  as 
described  in  the  baseline  problem. 

Given  the  continuous  time  plant 

_y _ ks  +  b 

a  corresponding  discrete-time  model  is  given  by 

y  kz\bz^ 

^ref  1  +  a^Z  ’  +  a^z  ^ 

The  system  can  then  be  written  in  a  recursive  fashion  as  a  difference 
equation 
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y(n)  +  a^yln-l )  +  a,y{n-2)  =  ky,g,(n-1 )  +  by^g,(n-2) 
y(n)  =  -a^ytn-l )  -  a,y(n-2)  +  ky^,(n-i )  +  by^g,(n-2) 
y(n)  =  (t>^0 
where 

=  y(n-1 ),  y(n-2),  y,e,(n-1 ).  y,e,{n-2) 
e^^f-a^,  -a^,k,  6] 

The  vector  4)  is  a  regression  vector  of  the  inputs  and  outputs  of  the  system. 

The  vector  6  is  a  parameter  vector  which  contains  the  system  parameters 
(coefficients  of  the  open-loop  system  transfer  function).  We  can  now  generate 
a  new  system  defined  by 


0(n+1)  =  0(n) 
y(n)  =  /(n)0(n) 


where  9,4)  are  exactly  as  before. 

The  appropriate  RLS  estimator  can  simply  oe  derived  as  a  Kalman  filter 
for  the  system  above.  The  corresponding  Kalman  filter  equations  are 


0(n)  =  0(n-1)  +  Kg{n)[y(n)  -(t)^(n)0(n-1)) 

Ke(n)  =  P(n)<tKn)  =  P(n-1  )(t>(n)tXI  +  (D^(n)P(n-l  )^n)] 

P(n)  =  [\-K^(nW(n)]P(n-m 


again,  where  4)  is  a  vector  of  inputs  and  outputs  and  0  is  the  vector  of 
parameters  which  are  to  be  estimated.  The  second  and  third  equations 
above  give  the  optimal  gain  Kg  of  this  estimator  in  terms  of  the  covariance 

matrix  P,  both  of  which  are  calculated  accordingly.  The  value  X  is  defined  as 
the  exponential  "forgetting  factor.”  This  forgetting  factor  is  used  to  discount 
older  data  which  becomes  invalid  if  the  parameters  vary  as  a  function  of  time. 
For  X  =1  no  data  is  discarded,  as  X  approaches  0  older  data  is  discounted 
more  quickly. 

The  use  of  a  discrete-time  plant  model  in  the  estimation  algorithm  brings 
along  with  it  a  subtle  but  very  important  point,  and  that  is  the  actual  values 
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being  estimated  are  the  coefficients  of  the  discrete-time  transfer  function  and 
not  the  continuous  time  transfer  function.  The  effect  of  transforming  a 
continuous  time  system  into  a  discrete  time  model  is  the  mathematical 
equivalent  of  mapping  the  entire  complex  plane  into  the  unit  circle,  as  is 
readily  recalled  from  elementary  digital  control  theory  (see  Figure  5.1).  The 
net  result  of  this  operation  appears  in  the  accuracy  of  the  estimates 
achievable  by  any  of  the  estimator  algorithms.  Since  such  a  very  large  area 
IS  compressed  into  the  unit  circle,  it  is  immediately  obvious  that  estimating  the 
parameters  to  within,  say  10  percent,  of  their  actual  values  in  the  discrete 
domain  certainly  does  not  mean  that  their  values  are  known  to  within  10 
percent  in  the  continuous  domain. 

As  a  further  example  of  this  problem  consider  the  general  continuous 
time  characteristic  equation 


(S  p,)(s  +  p^) . (s  +  pj 

which  has  a  corresponding  discrete-time  equivalent  (using  the  direct  pole¬ 
mapping  procedure  outlined  in  [10]) 


■P<T 

(z-e  )(Z' 


,  -ftT 

(Z-e  ) 


where  T  is  the  sampling  period.  Notice  that  as  T  approaches  0  (or  the  faster 

the  system  is  sampled),  that  e'P^  approaches  1.  and  the  discrete-time 
characteristic  equation  approaches 


(1  -z-b" 

which  results  in  all  of  the  roots  of  the  characteristic  equation  becoming 
concentrated  about  the  point  (-1 ,0)  in  the  discrete  domain.  This  means  that 
as  the  system  is  sampled  faster  that  the  discrete-time  characteristic  equation 
will  tend  towards 


(1  -z-"')'^ 

regardless  of  the  actual  underlying  continuous  time  characteristic  equation 
which  the  discrete  time  model  represents. 

The  practical  implication  of  this  result  is  that  if  very  fast  sample  rates  are 
to  be  used,  then  the  number  of  significant  digits  in  the  estimates  required  will 
be  quite  large.  As  an  example  consider  the  second  order  continuous  time 
characteristic  equation 
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s2  +  2s  +  1 


which  has  the  corresponding  discrete-time  model 

1  -  (e'^-r  e'^^)z'''  -t-  (e'^^)z'^ 


For  the  following  sample  periods  T,  this  equation  becomes 

1  -  0.5z  ’ +  0.05z‘^,  T=1.0 
1  -  1.72z'^  +  0.74z'^.  T=  0.1 
1  -  1.98Z  '  +  0.98z  ^,  T=  0.01 

The  net  result  is  that  as  the  system  is  sampled  faster,  the  discrete-time 
coefficients  become  less  distinguishable.  This  is  a  totally  counterintuitive 
result  and  is  one  that  is  not  at  all  apparent  to  designers  who  are  not  intimately 
familiar  with  discrete  recursive  parameter  estimation  or  STR  adaptive  control 
design. 

Fortunately  there  is  an  alternative  to  performing  estimation  on  discrete- 
time  models  of  continuous  systems,  and  that  is  to  estimate  the  continuous 
time  parameters  directly.  The  delta  operator  method  of  [9]  is  exactly  suited  for 
this  application  and  can  be  implemented  using  the  same  estimator  routines 
di.scussed  earlier.  Using  the  delta  operator  approach  it  is  true  that  the  faster 
the  system  is  sampled  the  more  accurate  the  estimates  become,  which  is  the 
desired  effect  of  fast  sampling.  In  fact,  it  is  not  unusual  to  be  able  to  estimate 
the  continuous  time  parameters  to  within  1  percent  of  their  nominal  values 
under  suitable  conditions. 

In  view  of  all  of  the  preceding  discussion,  the  primary  limitations  of  STR 
adaptive  control  can  be  summarized  by  considering  two  very  important  cases. 
The  first  limitation  is  the  need  for  persistency  of  excitation  in  order  to  perform 
any  kind  of  parameter  estimation.  Consider  again  the  system  equation 


y(n)  =  -a,^y(n-1 ) ...  -  a,y(n-m)  +  b,^y^^,(n-1 ) ...  +  b^y^^,(n-m) 


If  the  input  y^g,  =  0,  then  it  is  obvious  that  no  information  about  b^,  b2,...,  will 
be  available.  Also,  if  y^g,  =  0  and  the  system  is  at  rest  (y  =  0),  then  no 
information  concerning  a, ,  a2......  will  be  available  either.  This  results  in 

the  requirement  imposed  throughout  STR  design  literature  that  the  input  must 
be  "sufficiently  rich"  in  terms  of  excitation  before  STR  adaptive  control  can  be 
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successfully  employed. 


A  second  and  more  debilitating  limitation  of  STR  adaptive  control  is  the 
inability  of  any  STR  algorithm  to  cope  with  unmeasurable  external 
disturbances  which  are  similar  in  frequency  content  to  the  reference  input. 
This  can  be  easily  demonstrated  by  considering  the  first-order  example 
shown  in  Figure  5.2.  The  reference  input  is  u,  w  is  an  unmeasurable  external 
disturbance,  and  a  is  the  plant  parameter  which  has  a  nominal  value  a  =  1 
Take  for  example  the  case  when  u  =  1  and  the  output  y  is  observed  to  be  the 
value  y  =  2.  There  is  no  way  to  discern  from  this  limited  set  of  input  and 
output  data  whether  the  variation  in  y  from  it's  nominal  value  of  1  to  the  actual 
value  of  2,  is  due  to  a  perturbation  in  a  (i.e.  the  actual  value  of  a  might  be  0.5 
as  opposed  to  the  nominal  expected  value  of  1),  or  whether  in  fact  an  external 
disturbance  of  w  =  1  exists.  The  fundamental  problem  is  that  there  is  one 
equation  which  has  two  unknowns.  One  could  arbitrarily  postulate  that 
certain  observed  output  effects  are  caused  by  either  the  external  disturbance 
or  the  parameter  perturbation,  but  there  is  no  way  to  uniquely  determine  a 
cause  and  effect  relationship.  Most  of  the  current  STR  design  algorithms 
attempt  to  circumvent  this  problem  by  isolating  the  input  and  output  signals 
from  the  external  disturbances  in  frequency  by  band-pass  filtering  itiese 
signals  appropriately,  but  there  is  no  cure  for  the  case  when  the  external 
disturbance  has  a  similar  frequency  content  to  the  input  signal. 

As  an  example  consider  the  following  system 

y  _  ks  -t-  b 

V  2 

ref  s  +  a^s  +  a, 


where  the  nominal  values  of  a,  and  a2  are  1 .0  and  2.0  respectively.  Figure 
5.3  shows  the  estimates  of  the  parameters  when  no  disturbance  is  present  (w 
=  0),  and  Figure  5.4  shows  the  estimates  produced  by  the  exact  same 
algorithm  when  a  constant  external  disturbance  is  present  (w  =  1).  The 
estimator  used  in  this  case  was  a  RLS  with  exponential  forgetting  type.  It  is 
immediately  obvious  that  the  presence  of  an  external  disturbance  severely 
limits  the  performance  of  the  parameter  estimator. 

in  view  of  the  limitations  and  design  considerations  associated  with  STR 
adaptive  control  as  generally  outlined  in  this  report,  it  is  a  logical  conclusion 
that  this  particular  method  of  adaptive  control  is  not  at  all  suitable  for  the  type 
of  realistic  problems  embodied  by  the  benchmark  problem  considered  herein. 
This  in  no  way  suggests  that  STR  adaptive  control  does  not  have  merit  as  a 
viable  design  technique,  but  it  does  strongly  indicate  that  such  an  approach  is 
inappropriate  for  the  type  of  aerospace  problems  of  interest  in  this  study. 
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ENTIRE  LEFT-HALF 


Figure  5.1.  Mapping  of  Continuous  Domain  into  Discrete  Domain 


Figure  5.2.  Effects  of  External  Disturbances  on  Simple  Estimation 
Problems 
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Figure  5.3.  Continuous  Time  Parameter  Estimation,  No  External 
Disturbances 


Continuous  Time  Parameter  Estimation,  With  External 
Disturbances 


6.0  Model  Reference  Adaptive  Control 

The  second  major  area  of  adaptive  control  theory  is  known  as  Model 
Reference  Adaptive  Control  (MRAC).  Implementation  is  achieved  by 
constructing  an  ideal  model  which  exhibits  the  desired  response  to  some 
command  reference,  then  comparing  the  actual  plant  output  to  the  ideal 
model  output  to  obtain  an  error  signal.  The  error  signal  is  used  to  modulate 
the  feedback  gains  of  the  controller  so  as  to  eventually  drive  the  difference 
between  the  ideal  response  and  the  actual  response  to  zero  (see  Figure  6.1). 

The  origins  of  this  type  of  control  strategy  can  be  traced  back  to  the  "MIT 
rule"  developed  at  the  Massachusetts  Institute  of  Technology  in  the  late 
1950's.  The  MIT  rule  was  a  straightforward  approach  that  suggested  that  the 
controller  gains  should  be  changed  in  proportion  to  the  gradient  of  the  e^ror 
signal.  Very  simply  put,  the  larger  the  error  the  larger  the  gain  required,  and 
the  smaller  the  error  the  smaller  the  gain  required.  There  have  been  many 
imorovements  and  iterations  made  upon  the  basic  concept  over  the  last  three 
decades,  but  most  of  the  current  MRAC  techniques  still  contain  these  same 
basic  concepts. 

A  detailed  discussion  of  the  specific  MRAC  design  algorithm  used  m  this 
study  can  be  found  in  [1 1],  Sample  applications  are  discussed  in  [12]  and 
[13],  An  outline  of  the  controller  design  development  is  summarized  in  the 
following  discussion. 

Given  a  reference  model  system 

Vm  “ 


and  the  plant  equations 


=  ApXp  +  BpUp 


yp  ~ 


Define  the  error  as 


®  ~  ym '  yp 


A  reference  vector  is  constructed 
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where  e  is  the  error  as  defined  above,  is  the  ideal  model  state  vector,  and 
Urn  is  the  reference  input  (possibly  a  vector).  The  adaptive  control  iaw  is  given 
by 


u  =  kar 


where 


+  k| 


and 


kp  =  er^T=e(e.x^.  uJT 
k,  =er^T=e[e.Xn,.ujT 


and  the  gain  kg  can  be  rewritten  as 

where 

k^=e2T+ je^Tdt 

^=eu^T  +  Jeu^Tdt 
The  adaptive  control  law  can  then  be  written  as 


u  =  eke  +  k^Xe,  +  k^u^ 


As  a  specific  example  consider  the  baseline  problem  which  gives  the 
following  plant  model 
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ks  +  b 


_y 

y  2 

^ref  S  +  328  + 


This  plant  can  be  written  in  the  following  state-space  canonical  form 

^  ■  Vi  ^ 

x2  =  -aix,  +  bup 

y  =  x, 


The  reference  model  for  this  system  represents  the  desired  closed  loop 
response  in  transfer  function  form,  and  is  given  by 

y  _  0.5s  ^  1 
+  2S  +  1 


which  can  be  written  in  state-space  format  as 

X  .  =  -2x  +  X  ^  +  0.5u„ 


ml 


ml  m2 


X  ^  =  -  X  ,  + 
m2  m1  ^ 


ynr’^ml 


Now  define 


and  the  reference  vector 


e  =  Vm  -  y 


and 


where 


u  =  kar 


40 


where 


k  =er^T=e[e,  X  ,,x  ..  u  JT 

P  ^  ’  ml’  m2  fTV 


k,  =  er"T=e[e,x^^.x_^2’^-]T 


then 


i'  ^  (k  k  k  k  ] 

c*  t  “e»  *'x, »  ^ui 


kg  =  e^T  +  Je^T  dt 
k,  =  ex  T  +  lex  T dt 


ml 


h 


ml 


^  dt 


Finally,  the  complete  adaptive  control  law  is 


u  =  ek-  +  k,  X  +  k^  X  „  +  k^.u„ 
®  *1  ml  *2  rn2  ^  ^ 


A  block  diagram  of  the  complete  closed  loop  system  is  shown  m  Figure  6.2, 


T  and  T  are  weighting  matrices  to  be  chosen  by  the  designer.  The 
particular  choice  of  these  parameters  will  determine  the  adaptation  rate  and 
the  domain  of  adaptability  achievable  by  the  controller.  There  are  no 
analytical  procedures  available  to  aid  in  the  selection  of  these  weighting 
matrices,  even  though  they  are  critical  in  determining  the  operation  of  the 
closed  loop  controller.  The  most  straightforward  method  of  choosing  these 
weighting  matrices  is  through  an  empirical  study  of_observed  ideal  model 
following  performance  for  specific  choices  of  T  and  T. 

This  particular  adaptive  control  scheme  exhibits  very  good  performance 
characteristics  in  terms  of  maintaining  the  ideal  response  when  the  actual 
plant  is  subject  to  significant  uncertainties.  The  plots  shown  in  Figures  6.3, 
6.4,  6.5,  and  6.6  illustrate  the  perforrnance  capabilities  of  the  MRAC  controller 
as  a  function  of  the  choice  of  T  and  T.  In  each  case  the  reference  model  is 
described  by  the  reference  plant  above,  with  values  k  =  0.5,  b  =  1 .0,  a,  =  l  .0, 
a2  =  2.0.  The  actual  plant  has  parameters  k  =  0.5,  b  =  1 .0,  a^  =  1 .0,  and  a2  = 
0.0.  Notice  that  the  expected  (or  nominal)  plant  parameters  are  the  values 
given  by  the  ideal  model  description.  The  actual  value  of  a2  differs 
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Figure  6.2.  MRAC  Controller  Structure 


_ ~1 


i 


Figure  6.3.  MRAC  Controller  Performance,  T,T  =  C.01 

lP 

-1 


Figure  6.4.  MRAC  Controller  Performance,  T,T  =  0.1 
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Q  O'"  16  0  32  0  ^  48.0  64.0  ^.0 

O' 

Figure  6.5.  MRAC  Controller  Performance,  T,T  =  1.0 


o 

lP 


Figure  6.6.  MRAC  Controller  Performance,  T,T  =  10.0 
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significantly  from  the  expected  value.  In  this  case  32=  0.0  represents  a 
system  with  no  damping  even  though  the  controller  was  designed  expecting 
an  overdamped  response. 

It  IS  immediately  obvious  as  to  the  impact  that  the  choice  of  T  and  T  has 
upon  the  MRAC  performance  capabilities.  The  Figures  of  6.3,  6.4,  6.5,  and 
6.6  illustrate  the  effect  of  increasing  these  tuning  matrices  until  the  actual 
plant  output  is  made  to  coincide  exactly  with  the  ideal  response. 
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7.0  Classical  Controller  Design  Comparison 

In  addition  to  comparing  the  adaptive  control  designs  to  one  another,  it  is 
also  important  to  compare  the  designs  with  existing  accepted  technology. 

This  is  really  the  only  way  to  judge  the  relative  merits  and  limitations  of  these 
advanced  techniques.  In  this  particular  case  a  traditional  Proportional- 
Integral-Derivative  (PID)  compensator  was  chosen  for  comparison  purposes. 

A  PID  controller  was  placed  in  series  with  the  nominal  plant  nf  the 
benchmark  problem  and  was  tuned  to  give  the  desired  closed  loop  response. 
Figure  7.1  shows  both  the  open  loop  response  (no  PID  in  the  loop)  and  the 
closed  loop  response  (with  PID  in  the  loop)  for  the  nominal  plant. 

A  common  ^orm  of  the  PID  control  law  is 

k 

u  =  k  +  -^  +  k  s 

PS  d 

In  this  problem  the  phase  advance  supplied  by  the  derivative  term  was  not 
necessary  and  was  even  undesiiable,  therefore  the  gain  was  set  to  zero. 
This  in  effect  makes  the  control  law  proportional-plus-integral  (PI).  The  gains 
k,  and  kp  were  subsequently  chosen  to  provide  the  closed  loop  response 
shown  in  Figure  7.1.  The  resulting  gains  values  were 

kp=  I.O 
k.  =  0.55 

I 

k^  =  0.0 

a 
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0.30  0.70 


8.0  Performance  Results 

The  controller  designs  discussed  in  sections  4.0  through  7.0  were 
compared  in  terms  of  their  performance  robustness  capabilities  when 
subjected  to  significant  uncertainties  due  to  parameter  perturbations,  external 
disturbances,  and  unmodeled  dynamics.  The  dynamical  model  used  in  each 
of  the  designs  was  the  plant  described  in  the  benchmark  problem  of  section 
2.2. 

In  each  case  the  controllers  were  judged  in  terms  of  their  ability  to 
maintain  the  same  desired  ideal  plant  response  that  was  given  by  an  ideal 
model  of  the  closed  loop  system.  The  plant  dynamics  are  governed  by 


'l_  ks  +  b 
u  ~  2 

S  +  32^  +  3.^ 

The  ideal  closed  loop  response  is  given  by  the  plant  transfer  function  above 
with  the  following  set  of  nominal  parameters 


^ideal  _  0.5s  h-  1 
^ideal  +  2S  +  1 


Figure  8.1  shows  the  ideal  input  and  desired  ideal  output  for  the  system. 

A  performance  index  was  constructed  to  quantitatively  measure  the 
ability  of  the  controllers  to  maintain  the  ideal  response  when  subjected  to 
uncertainties  and  off  nominal  plant  parameters.  The  performance  index  was 
chosen  to  reflect  the  type  of  requirements  illustrated  in  Figure  3.1  of  section 
3.0.  In  order  to  measure  both  the  transient  behavior  and  the  steady-state 
tracking  capabilities  of  the  different  controllers,  the  following  performance 
index  was  used 


J  =  (j e^dt)  w  (lei) 
0 


where 


®  =  y,deai 


-y 


actual 


48 


This  performance  index  is  the  union  of  Integral-squared-error  (ISE)  and 
absolute  error  (AE),  where  the  AE  term  penalizes  the  transient  behavior  in 
terms  of  peak  deviation  away  from  the  ideal  response  and  the  ISE  term 
penalizes  poor  steady-state  tracking  performance. 

t 

8.1  Parameter  Perturbations 

The  ADAC  controller,  conventional  DAC  (polynomial  spline),  MRAC 
controller,  and  PID  controller  were  evaluated  in  terms  of  their  ability  to 
maintain  the  ideal  model  response  when  subject  to  parameter  perturbations 
away  from  the  expected  nominal  values.  Each  of  the  controller  designs  were 
outlined  in  sections  4.1 , 4.2,  6.0,  and  7.0  respectively.  Appendix  A  contains 
listings  of  the  simulation  programs  used  to  perform  these  comparisons.  In  the 
listings  are  the  detailed  equations  for  each  of  the  control  laws  used.  These 
simulations  are  written  in  ACSL  (Advanced  Continuous  Simulation 
Language)  and  the  equations  are  self  documenting.  Figure  8  2  shows  a 
block  diagram  of  the  test  setup  used  to  evaluate  the  controller  designs. 

In  each  case  the  controller  was  placed  in  the  closed  loop  with  the 
nominal  plant.  The  plant  parameters  were  then  varied  away  from  their 
nominal  values  and  the  performance  index  was  monitored.  The  parameters 
were  individually  varied  until  the  performance  index  exceeded  a  prescribed 
threshold.  By  individually  varying  the  parameters  the  entire  parameter  space 
was  explored  and  certain  areas  of  suitable  performance  were  established. 

The  threshold  on  the  ISE  term  of  the  performance  index  was  1 .0  and  the 
threshold  for  the  AE  term  was  0.2.  As  long  as  these  thresholds  were  not 
exceeded  the  controller  was  able  to  maintain  the  desired  closed  loop 
response  satisfactorily. 

Figures  8.3,  8.4,  and  8.5  show  plots  of  the  ADAC  controller,  Polynomial 
Spline  DAC,  and  PID  controllers  in  the  a^  -  83  parameter  space  for  the  fixed 
values  of  b  =  1 .0,  b  =  3.0,  and  b  =  1 0.0  respectively.  Note,  only  a^ ,  05,  and  b 
were  varied  in  this  study,  so  the  parameter  space  can  be  visualized  as  a  three 
dimensional  Euclidian  space.  The  effect  of  varying  b  was  to  change  the 
location  of  the  plant  zero,  which  indirectly  effects  the  gain  term  as  well, 
therefore  to  simplify  the  results  the  gain  parameter  k  was  maintained  at  the 
nominal  value  of  0.5  always.  Keep  in  mind  that  the  enclosed  regions  in 
Figures  8.3,  8.4,  and  8.5  depict  areas  of  unexpected  off-nominal  parameter 
variations  where  the  ideal  model  response  was  maintained.  This  distinction 
is  important  in  that  these  regions  indicate  the  performance  robustness  of  the 
various  controller  designs,  as  opposed  to  the  more  popular  but  less  important 
stability  robustness  which  is  often  cited  in  controls  literature. 

Examination  of  Figures  8.3,  8.4,  and  8.5  reveals  that  both  of  the  DAC 
controllers  are  far  superior  to  the  PID  controller,  as  was  expected.  Variations 
in  terms  of  orders  of  magnitudes  were  tolerable  for  the  ADAC  and  Polynomial 
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Figure  8.2.  Controfier  Test  Block  Diagram 


Figure  8.3.  Parameter  Space  Performance  Plots,  b=l 


5 


Figure  8.4,  Parameter  Space  Performance  Plots,  b=3.0 
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Figure  8.5.  Parameter  Space  Performance  Plots,  b=l0. 
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spline  DAC  controllers.  Notice  in  Figures  8.3  and  8.4  that  the  polynomial 
spline  DAC  seems  to  be  relatively  insensitive  to  variations  in  a^,  while  the 
ADAC  design  reflects  a  nice  symmetry  about  the  nominal  plant  parameters. 
This  observation  is  even  more  apparent  in  Figures  8.6  and  8.7  \^  lich  show 
the  ADAC  and  Polynomial  Spline  DAC  controllers  separately,  for  the  different 
values  of  b.  However,  in  both  Figure  8.5  and  8.7  it  is  evident  that  the 
Polynomial  Spline  DAC  is  very  sensitive  to  variations  in  b,  and  for  b  =  10.0 
the  region  of  acceptable  performance  has  collapsed  to  a  single  point  about 
the  nominal  plant  values.  The  ADAC  design  does  not  reflect  such  a 
sensitivity  to  b  variations,  and  Figure  8.6  indicates  that  even  for  the  case  b  = 
10.0  the  area  of  acceptable  a^  -  a2  variations  is  still  quite  large.  This  sei  of 
plots  indicates  that  if  the  regions  of  acceptable  performance  were  stacked  on 
top  of  each  other  in  a  three  dimensional  (a,,  a2,  b)  space,  that  the  resulting 
volume  of  acceptable  performance  for  the  ADAC  controller  would  be  larger 
than  the  corresponding  volume  for  the  Polynomial  Spline  DAC  controller.  If 
only  the  plots  of  Figures  8.3  and  8.4  were  considered,  it  would  appear  that  the 
polynomial  spline  DAC  controller  was  better  than  the  ADAC  controller,  but  if 
the  entire  volume  of  parameter  variations  is  considered  the  ADAC  controller 
is  better.  The  decreased  insensitivity  of  the  polynomial  spline  DAC  compared 
to  the  ADAC  for  certain  parameter  variations  shown  in  Figures  8.3  and  8.4  is 
due  to  the  increased  ability  of  the  polynomial  spline  model  to  estimate 
particular  characteristics  of  the  system  uncertainty  as  compared  to  the 
exponential  basis  functions  used  in  the  ADAC  controller.  Figure  8.8  shows 
the  PID  controller  performance  for  the  same  set  of  b  variations. 

The  results  of  Figures  8.6,  8.7,  and  8.8  are  somewhat  misleading  as  they 
indicate  the  controllers  better  performance  for  off  nominal  values  of  b.  This  is 
an  artifact  of  the  choice  of  performance  index.  For  the  value  of  b  =  3.0  the 
performance  region  is  larger  than  that  for  b  =  1 .0  or  b  =  10.0  since  neither  the 
absolute  error  (AE)  or  the  integral-squared-error  (ISE)  terms  in  the 
performance  index  dominates.  Rather,  both  terms  are  relatively  equal  in 
magnitude.  For  the  cases  of  b  =  1 .0  and  b  =  10.0,  either  the  ISE  or  AE  term 
dominates  the  performance  index  and  the  boundary  is  defined  by  primarily 
one  term  or  the  other.  The  particular  combination  of  parameter  variations  for 
the  case  b  =  3.0  results  in  the  boundary  being  defined  by  a  combination  of 
both  terms  equally.  Another  choice  of  performance  index  (strictly  ISE  alone 
or  AE  alone  would  not  show  the  performance  improving  as  b  moved  away 
from  nominal).  The  performance  index  chosen  is  still  valid  in  terms  of  making 
comparisons,  however,  since  all  controllers  were  judged  using  the  same 
criteria. 

The  dynamical  order  of  the  ADAC  and  Polynomial  Spline  DAC 
controllers  was  the  same  (6th  order).  This  is  in  contrast  to  the  PID  controller 
which  has  a  dynamical  order  of  one.  The  poles  of  the  composite  observers  for 
both  the  ADAC  and  Polynomial  Spline  designs  were  identically  located  at 
(-3.0)  in  the  complex  plane  (all  poles). 
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Figure  8.9  shows  a  similar  set  of  plots  for  the  MRAC  controller  design. 
The  performance  area  are  plotted  for  several  values  of  the  tuning  matrices  T 
and  T,  It  is  obvious  that  these  tuning  parameters  could  be  increased  until  the 
regions  of  acceptable  performance  covered  practically  the  entire  parameter 
space.  For  this  reason  it  is  very  difficult  to  compare  the  the  results  of  the  DAC 
controllers  with  the  MRAC,  since  wejiave  no  way  of  gauging  a  set  of  realistic 
values  for  the  tuning  matrices  T  and  T. 

A  set  of  typical  transient  responses  and  the  associated  control  signals  for 
the  ADAC,  MRAC,  and  PID  controllers  are  shown  in  Figures  8.10,  8.1 1  and 
8.12  respectively.  These  responses  are  indicative  of  the  types  of  responses 
existing  when  the  system  is  very  close  to  one  of  the  boundaries  shown  in  the 
parameter  space  plots  above.  In  other  words,  these  transient  responses 
represent  "worst-case"  acceptable  responses.  For  the  ADAC  controller  tne 
particular  parameter  values  were  b  =  1 .0,  a-i  =  -5.0,  ag  =  -5.0.  For  the  MRAC 

controller  these  parameters  were  b=1 .0,  a^  =  1 0.0,  32  =  -5.0,  T  =  1 0  0.  For 
the  PID  controller  these  parameters  were  b=1 .0,  a-]  =  i  .0,  32  =  -0.5. 

8.2  External  Disturbances 

The  baseline  closed  loop  system  was  subjected  to  external  disturbances 
as  shown  in  Figure  8.13. 

Each  of  the  adaptive  controller  designs  ADAC,  MRAC,  Polynomial 
Spline  DAC,  and  the  classical  PID  design  use  some  type  of  integral  feedback 
in  their  control  laws.  This  means  that  each  of  these  controllers  will  be 
capable  of  some  d^gr^o  of  disturbance  •■ejection.  The  actual  degree  of 
rejection  possible  is  a  function  of  the  controller's  response  lime  (or  inversely, 
the  controller  bandwidth). 

The  response  time  of  the  DAC  controllers  is  set  by  the  choice  of  the 
observer  poles.  The  deeper  the  poles  are  placed  into  the  ieft-half  of  the 
complex  plane,  the  faster  the  response,  and  the  greater  the  bandwidth.  For 
the  MF^AC  controller,  the  response  time  is  set  by  the  choice  of  tuning  matrices 
T  and  T.  The  greater  the  value  of  T  and  T,  the  faster  the  response,  and  the 
greater  the  bandwidth.  For  the  PID  controller  the  disturbance  response  is  set 
by  the  choice  of  the  integral  gain  k,. 

To  demonstrate  the  disturbance  rejection  capabilities  of  each  controller, 
they  were  individually  subjected  to  a  set  of  disturbances  that  covered  the 
entire  possible  spectral  space  of  the  controller's  bandwidth.  This  disturbance 
set  was  comprised  of  two  basic  signals.  The  first  disturbance  signal  was  a 
single  sinusoid,  varied  in  frequency  from  dc  to  beyond  the  controller 
bandwidth.  This  single  spectral  line  swept  through  the  entire  controller 
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Figure  8.12.  PIO  Controller  Transient  Response  and  Control 
Signal 
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Figure  8.13.  External  Disturbance  Inputs  to  the  System 
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spectrum  represents  every  single  disturbance  frequency  possible.  The 
second  disturbance  signal  was  wideband  white  noise.  The  spectral  content 
of  white  noise  is  a  continuum  of  frequencies  which  cover  the  controller 
bandwidth.  This  choice  of  disturbance  represents  every  possible 
combination  of  dioturbance  frequencies.  These  spectral  characteristics  are 
shown  graphically  in  Figure  8.14.  Time  domain  examples  of  these 
disturbance  signals  are  shown  in  Figures  8.15  and  8.16. 

The  response  of  the  ADAC,  MRAC,  and  PID  controllers  to  these 
disturbances  are  shown  in  Figures  8.17,  8.18,  and  8.19  respectively.  Figure 
8.17  illustrates  the  ability  of  the  ADAC  controller  to  maintain  the  ideal 
response  characteristics  in  the  presence  of  any  disturbance  which  has  a 
spectral  content  that  is  within  the  controller  bandwidth.  The  disturbance 
rejection  capabilities  of  the  ADAC  controller  are  not  surprising  since  the  entire 
concept  of  DAC  control  was  originated  for  this  particular  purpose.  The 
disturbance  rejection  performance  can  be  better  appreciated  by  comparing 
Figure  8.17  with  Figure  8.19  which  shows  the  PID  controller  performance. 

The  PID  has  limited  capability  due  to  the  single  integrator  available  in  the 
controller.  In  fact  it  can  be  shown  analytically  that  the  only  type  of  disturbance 
which  can  be  completely  canceled  by  this  PID  controller  is  a  constant 
disturbance 

The  actual  outputs  (control  signals)  of  the  ADAC,  MRAC,  and  PID 
controllers  are  shown  in  Figures  8.20,  8.21,  and  8.22  for  comparison 
purposes.  These  are  the  control  signals  for  the  case  of  the  single  sinusoidal 
disturbance. 

8.3  Unmodeled  Dynamics 

Unmodeled  dynamics  were  introduced  into  the  baseline  system.  The 
unmodeled  dynamics  can  be  separated  into  two  categories:  (1)  multiplicative 
and  (2)  additive.  The  distinction  between  these  two  types  of  unmodeled 
dynamics  has  a  significant  impact  on  the  control  system  design  and 
performance. 

Consider  the  case  of  multiplicative  unmodeled  dynamics  shown  in 
Figure  8.23.  This  Figure  shows  the  characteristic  pole  locations  in  the 
complex  frequency  plane  for  a  system  with  eight  total  poles,  two  of  which  are 
modeled  in  the  assumed  dynamics.  Those  poles  lying  well  to  the  left  of  the 
modeled  poles  are  assumed  to  be  "insignificant"  in  terms  of  their  effect  on  the 
overall  system.  The  fact  that  these  poles  do  not  contribute  significantly  to  the 
response  of  the  system  allows  them  often  to  be  disregarded  during  the  design 
process.  This  is  standard  practice  in  control  system  design  and  is  a  well- 
accepted  and  valid  approach.  There  are  four  other  unmodeled  poles  in  this 
example  which  lie  to  the  right  of  the  modeled  poles.  The  fact  that  they  lie  to 
the  right  in  the  complex  plane  means  that  these  poles  will  contribute 
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Figure  8.14.  Spectral  Characteristics  of  External  Disturbances 
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Figure  8.20.  ADAC  Control  Signal,  Sinusoidal  Disturbance  Input 
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Figure  8.21.  MR  AC  Control  Signal,  Sinusoidal  Disturbance  Input 
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Figure  8.22.  PID  Control  Signal, 
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Figure  8.23.  Characteristic  Poles  of  Muitipllcitive  Unmodeled 
Dynamics 
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Significantly  to  the  overall  response  and,  in  fact,  will  dominate  the  system 
response.  Designing  a  control  system  based  upon  particular  modeled 
dynamics  while  neglecting  dominant  unmodeled  poles  is  tantamount  to 
designing  a  controller  based  upon  an  invalid  reduced-order  model  (ROM). 
The  proolems  associated  with  control  system  design  based  upon  ROM's  is 
well  known  and  documented  [14], 

None  of  the  adaptive  control  designs  (Polynomial  DAC,  ADAC,  and 
MRAC)  were  able  to  maintain  the  performance  requirements  in  the  presence 
of  dominant  multiplicative  unmodeled  dynamics.  This  really  is  a  somewhat 
impractical  problem  in  that  it  is  desired  to  design  a  control  system  for  a  plant 
in  which  the  dominant  characteristics  are  unknown,  and  an  invalid  ROM  is 
used  instead.  Not  even  adaptive  control  offers  a  solution  to  this  type  of 
problem. 

The  case  of  additive  unmodeled  dynamics  is  a  completely  different 
situation.  Additive  unmodeled  dynamics  act  as  effective  external 
disturbances  on  the  closed  loop  system.  It  is  apparent  that  the  controller's 
performance  should  be  similar  to  that  outlined  in  section  8  2.  This  is  in  fact 
the  case  and  Figures  8.24  and  8.25  show  an  example  of  this  performance  for 
the  ADAC  controller,  when  a  pair  of  undamped  unmodeled  poles  are  added 
to  the  closed  loop  system  as  shown  in  Figure  8.26  (the  frequency  of  these 
undamped  poles  was  0.25  cycles  per  second).  Figure  8.24  illustrates  the 
response  in  the  open  loop  case  (no  ADAC  control),  and  Figure  8.25  shows 
the  closed  loop  performance  (with  ADAC  control).  The  ADAC  system 
performs  very  well  in  the  presence  of  this  particular  additive  unmodeled 
dynamical  uncertainty. 
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Figure  8.26.  Transfer  Function  Format  of  Onmodeied 
Dynamics 
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9.0  Summary 

Disturbance  Accommodating  Control  principles  were  applied  to  the 
design  of  a  baseline  control  system.  The  resulting  DAC  controller  designs 
were  evaluated  in  terms  of  their  ability  to  maintain  an  ideal  model  response 
when  the  closed  loop  system  was  subjected  to  significant  uncertainties. 

These  uncertainties  included  parameter  perturbations,  external  disturbances, 
and  unmodeted  dynamics.  The  DAC  designs  were  compared  to  other 
adaptive  techniques  including  Self-Tuning  Regulator  and  Model  Reference 
Adaptive  Control  concepts.  These  designs  were  also  compared  to  a  classical 
Proportional-Integral-Derivative  controller.  All  of  the  adaptive  controller 
designs  were  judged  relative  to  each  other  and  relative  to  the  PID  design,  in 
terms  of  performance,  design  complexity,  and  confidence  of  operation 
(reliability). 

For  the  case  of  parameter  perturbations,  performance  results  were 
calculated  and  used  to  determine  specific  regions  in  a  three  dimensional 
parameter  space  where  the  controllers  were  able  to  maintain  ideal  response 
characteristics.  These  regions  were  centered  about  some  nominal  set  of 
parameter  values,  and  their  size  indicates  the  amount  of  parameter 
uncertainty  allowable  for  the  various  controller  designs. 

For  external  disturbances  it  was  shown  that  the  disturbance  rejection 
capability  was  a  function  of  the  bandwidth  of  the  controller  and  the  spectral 
content  of  the  disturbance.  This  is  the  case  for  each  of  the  controller  designs 
which  use  some  type  of  integral  feedback.  Disturbance  rejection  was  verified 
by  subjecting  the  closed  loop  system  to  a  class  of  external  disturbances 
which  spanned  the  complete  spectral  space  possible.  This  set  of  external 
disturbances  was  composed  of  two  basic  components:  a  single  sinusoid 
varied  in  frequency  to  represent  every  single  possible  disturbance  frequency, 
and  broad  band  white  noise  to  represent  every  possible  combination  of 
disturbance  frequencies. 

Unmodeled  dynamics  were  classified  into  two  categories;  multiplicative 
and  additive.  Additive  unmodeled  dynamics  were  shown  to  have  much  the 
same  effect  as  external  disturbances  and  compensated  for  much  in  the  same 
manner  by  the  various  controllers.  Multiplicative  unmodeled  dynamics  were 
shown  to  present  a  problem  to  all  of  the  adaptive  controller  designs,  in  that 
none  of  the  controllers  were  able  to  maintain  performance  margins  in  the 
presence  of  dominant  multiplicative  unmodeled  poles. 

The  ADAC  controller  designs  were  shown  to  posses  very  good 
performance  robustness  in  the  presence  of  each  of  the  uncertainties 
discussed.  In  terms  of  parameter  perturbations  the  ADAC  controller  was  able 
to  withstand  variations  of  one  order  of  magnitude  for  both  the  numerator  and 
denominator  coefficients  of  the  nominal  plant  given  in  the  baseline  problem. 
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This  represents  a  vast  increase  in  capability  over  the  classical  PID  design, 
which  was  expected.  It  should  be  reiterated  that  all  of  the  results  shown  are 
in  terms  of  performance,  which  was  defined  as  the  ability  of  the  controller  to 
maintain  an  ideal  model  transient  and  steady-state  tracking  response.  This 
distinction  further  magnifies  the  regions  of  allowable  parameter  uncertainty 
shown  in  the  Figures  of  section  8.1. 

When  design  complexity  and  confidence  of  operation  are  considered  as 
crueria  along  with  performance  robustness,  then  it  becomes  very  obvious  that 
the  DAC  based  designs  are  superior  to  the  other  adaptive  control  techniques. 
The  MRAC  system  exhibited  very  good  performance  characteristics  for  certain 
choices  of  tuning  parameters,  but  the  inherent  nonlinearity  of  the  resulting 
controller  tends  to  limit  the  applicability  of  this  method.  The  only  tools 
available  to  analyze  and  measure  the  reliability  and  robustness  of  the  closed 
loop  system  are  repeated  simulation  and  Liapanov  stability  analysis.  Tiie 
resulting  DAC  based  controllers  are  linear  time-invariant  systems,  which 
means  that  all  of  the  standard  linear  analysis  tools  (Nyquist,  Bode,  root  locus, 
etc.)  are  siiii  applicable.  Again,  the  STR  approach  resulted  in  a  nonlinear 
closed  loop  system  which  was  subject  to  the  same  analysis  problems  as  the 
MRAC.  Even  beyond  the  issue  of  nonlinearity,  however  is  the  fundamental 
limitation  of  the  estimation  of  system  parameters  in  the  presence  of  external 
disturbances.  Figure  9.“  shows  a  matrix  which  details  these  comparisons. 
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Figure  9.1.  Relative  Comoarisoii  of  Adaptive  Control  Design 
Techniques 
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APPENDIX 


Computer  listings  of  controller  evaluation  simulations  written  in  ACSL 
code. 


A.1  Adaptive  DAC  (ADAC)  Simulation 


PROGRAM 

INITIAL 

LOGICAL  LDAC.LPULSE.LAOAC.LTVUF 
CONSTANT  A1IC-1.A2IC-2.AN1-1 .AN2-2 
CONSTANT  KIC-.5.BIC-1 .KN-.S.BN-l 
CONSTANT  MAl-0.  .MA2-0.  .MB-a. ,MK-0. 
constant  <01-14. 21.  <02-75.33 
CONSTANT  <03-351.5.  <04--281, 

CONSTANT  <05-18..  <06--782. 

CONSTANT  XlIC-0..  X2IC-0. 

CONSTANT  XMlIC-0..  XM2IC-0. 

CONSTANT  CHIIC-0..  E2HIC-0. 

CONSTANT  ZAlIC-0..  2A2IC-0. 

CONSTANT  ZlIC-0..  Z2IC-0. 

CONSTANT  PULSEP-20. 

CONSTANT  TSTOP  -79 . 99 . I SEMaX- 1 . . EMAX- . 5 
CONSTANT  WSINE-0. .WPULSE-0. .WNOISt-0. .WCONST-0. 
CONSTANT  MWPREQ-0. 1 .WPULSP-40. .WPULSW-20. 
CONSTANT  WMEAN-0. .WSTOEV-1 . .NOISEW-3. 

CONSTANT  LTVWF -. FALSE . 

CONSTANT  LOAC-.TRUE. . LPULSE- . TRUE .. LADAC- . TRUE . 
PULSEW-PULSEP/2. 

U  -  0. 
ch;-0. 


OVNAMIC 
OER : VATI VE 

CIN’ERVAL  CINT  -  0.05 
MAX'ERVAL  MAXT  -  0.01 
NSTE°S  -  1 


INPUT 


CELTA  -  RSV! LPULSE.PULSE(0. .PULSEP.PULSEW) , 1 . ) 

PARPMETER  VARIA'ICNS 

A1  -  MAl-T  ♦  AlIC 
A2  *  MA2*’"  *  A2IC 
3  -  MB'T  *  BIC 

<  -  MK'T  *  <:c 

■  CIS'UPBANCE  MODEL 

U  -  ws  INE-S  :  N  (  RSW(  LTVWF  .  Z«3  .  U15*MUFREQ-T  .  2-3  1  4  :  S  • '4UF  R  £0  )  -  T  )  ♦  ... 

WOU^SE*PUlSE(0. .WPULSP .WPULSW)  *  ... 

WNO:SE*OU( { 1/NOISEW) .WMEAN.WSTDEV)  *  ... 

WCONST 

SLANT  MODEL 

X 1 D  •  X2 

X  ;  .  ;  N  '  E  G  t  X  D  .  X  1 1  c  ) 

x:20  •  -Al*x;  -  A2-X2  *  u  *  U 
X2  -  1N''EG(X2D.X2IC) 

V  -  x-':b/<)*xi  ♦  x2) 

OBSERVEP  ■  ,  , 

A-  1 


e  .  DELTA  -  Y 
INOV  -  -  y  -  CHI 

CHID  -  E2H  ♦  <01*IN0V 
CHI  •  INTEG(CHID.CHIIC) 

E2H0  -  -ANl»CHI  -  AN2*E2H  -  AN1«DELTA  ♦  2A1H  -  ZIH... 

-BN"(UA  +  UD)  ♦  K02*IN0V 
E2H  -  INTEG(E2HD,E2HIC) 

2A1HD  •  ZA2H  +  K03*IN0V 
ZAIH  -  INTEG(ZAIHD.ZAIIC) 

ZA2H0  •  -ANl«ZAiH  -  AN2"ZA2H  ♦  K04«IN0V 
ZA2H  -  INTEG(ZA2H0.ZA2IC ) 

ZlHD  ■  Z2H  ♦  KOS'INOV 
ZIH  •  INTEGIZIHD.ZI IC) 

Z2HD  -  K06»INOV 

Z2H  -  INTEG(Z2HD,Z2IC) 

YH  •  -CHI 

•  CONTROL 

UP  -  (AN1/BN)*0ELTA 

UD  -  -  RSW( LDAC . { 1/BN )*Z1H ,0. ) 

UA  ■  RSW(LADAC. ( l/8N)•2AlH,^. ) 

U  -  (UP  ♦  UA  ♦  UD) 

•  IDEAL  MODEL  " 

XMID  -  XM2 

XMl  »  INTEGIXMID.XMI IC) 

XM20  -  -ANl'XMl  -  AN2*XM2  DEL^A 
XM2  •  I NTEG ( XM2D .XM2 IC ) 

YM  •  KN* ( { BN/KN  ) ‘XMl  +  XM2 ) 

■  INTEGRAL  SQUARED  ERROR  * 

ERROR  •  YM  -  Y 

ISE  •  INT£G(ERROR*'2.0. ) 

•  TERMINATION 

TERMT  ((ISE  .GE.  ISEMAX)  ,OR.  (ABSf'''’  .CE.  EMAX  > 

•OR.  <T  .GE.  TSTOP)) 

END  S“OF  DERIVATIVE" 


END  S’OF  DYNAMIC 


A.X  Polynomial  DAC  (ADAC)  Simulation 


PROCRAM 

INITIAL 

LOGICAL  LDAC.LPULSE 
CONSTANT  A1 IC-1 .A2IC-2.AN1-1  ,AN2«2 
CONSTANT  KIC-.5.BIC-1 .KN-.S.BN-l 
CONSTANT  MAl-0. .MA2-F. 

CONSTANT  KEl-1.846.  KE2-28.73 
CONSTANT  <E3-211.32,  <E4-695. 

CONSTANT  K:E5-n54..  KEB-TBl. 

CONSTANT  XlIC-0.,  X2IC-«. 

CONSTANT  XMlIC-iB..  XM2IC-iff. 

CONSTANT  XIHIC-B..  X2HIC-«. 

CONSTANT  ZIHIC-B..  22HlC-«..  23HIC-0..  Z4H1C-0. 
CONSTANT  PULSEP-2«. 

CONSTANT  TSTOP  •iggBB. 

CONSTANT  W-0. .LOAC-.TRUE. . L P ULSE- . TRUE . 
PULSEW-PULSEP/2 . 

U  •  0. 

END  S 

DYNAMIC 

DERIVATIVE 

CINTERVAL  CINT  •  0.05 
MAXTERVAL  MAXT  •  0.01 
NSTEPS  -  1 


INPUT  * 


DELTA  -  RSW<LPULSE,PULSE(0. .PULSEP.PULSEW) . 1 . > 

PARAMETER  VARIATIONS  " 

A1  -  MAl'T  +  AlIC 
A2  -  MA2*T  ♦  AZIC 
B  ■  MB*T  ♦  8IC 

V  •  MK*T  ♦  XIC 

PLANT  MODEL 

X  1  D  -  X2 

XI  ■  INTEGIXID.Xl IC) 

X2D  -  -A2«X2  -  Al*Xl  ♦  U  ♦  W 
X2  •  INTEG{X20.X2IC  ) 

Y  •  B»X1  +  K«X2 
OBSERVER 

DUM  -  Y  -  BN'XIH  -  KN*X2H 

XlHO  •  X2H  ♦  KEl’DUM 
XIH  ■  INTEG(XIHD.XIHIC) 

X2HD  •  -ANI'XIH  -  AN2*X2H  *  ZIH  ♦  U  ♦  RE2-DUM 
X2H  •  INTEG(X2HD .X2HIC) 

ZIHD  ■  22H  ♦  KED'-DUM 
ZIH  •  INTEG(ZIHD.ZIHIC) 


Z2HD 


Z3H  ♦  KE4*D'JM 


A-3 


-  INTEC(22HD.Z2HIC) 

23HD  -  ZiH  ♦  KE5*DUM 
23H  -  INTEG(Z3HD.Z3HIC) 

Z4hD  -  KE6«DUM 

2*H  •  INTEG(Z4H0 .Z4H1C ) 

•  CONT-  L 

U  ■  RSW(LDAC,-ZIH*(1/BN>«0ELTA,{I/BN)*0ELTA> 

•  IDEAL  MODEL  * 

XMID  -  XM2 

XMl  -  INTECCXMID.XMI IC) 

XM2D  -  -ANl-XMl  -  AN2*XM2  ♦  DELTA 
XM2  •  INTEC{XM2D,XM2IC) 

YM  -  KN*( (BN/KN )«XM1  *  XM2) 

•  INTEGRAL  SQUARED  ERROR  ■ 

ERROR  •  YM  -  Y 

ISE  ■  INTEG<ERROR**2.jB.  > 

•  TERMINATION 

TERMT  (USE  .GE.  ISEMAX)  .OR.  (ABS(ERROR)  .6E.  EMAX ) 
■OR.  IT  .GE.  TSTOP)) 

END  $*OF  DERIVATIVE* 

END  $*OF  DYNAMIC* 

END 


A-4 


A. 3  MRAC  Simulation 


PROGRAM 

iMiriAi 

LOGICAL  LADAPT.LPULSE.LERROR.LXMOO.LUM.LERRCH.LTWF .LTVWNS 
CONSTANT  A1 I C- 1 . A2 IC “2 . AN  1 " 1 . AN2-2 
CONSTANT  KIC-.5.B1C-1 ,KN«.S,BN-l 
CONSTANT  MAl-0. .MA2-0. .MB-*. .MK-*. 

CONSTANT  XlIC-*.,  X2IC-*. 

CONSTANT  XIMIC-*.,  X2MIC*-*. 

CONSTANT  CUC-*.  .C2IC-*.  , C3 1 C-* .  . C4 1C-*. 

CONSTANT  TN-20. .TNB-2*. 

CONSTANT  PULSEP-2*. 

CONSTANT  TSTOP  -79,99.  ISEMAX-1,.  EMAX-.l 
CONSTANT  WSINE-*. .WPULSE-0. .WNOISE-0. .WCONST-*. 

CONSTANT  MWrREQ-*. I .VPULSP-4*. .WPyLSW-2*. 

CONSTANT  WMEAN-ar. . WSTDEV- 1 . . NOISEW-3 • .MNOISE • .05 
CONSTANT  LTVWF-. false. . LTVWNS- . FALSE . 

CONSTANT  LADAPT-.TRUE. , L PULSE- . TRUE .. LERRCH- . TRUE . 

CONSTANT  L ERROR -.TRUE, . LXMOO- . TRUE .. LUM- . TRUE . 
PULSEU-PULSEP/2  . 

U  -  0. 

END  S 

DYNAMIC 

DERIVATIVE 

CINTERVAL  CINT  •  0. I 
MAXTERVAL  MAXT  -  0,01 
NSTEPS  ■  1 


■  INPUT  * 


DELTA  •  RSW( LPULSE .PULSE(0. .PULSEP .PULSEW) . 1 . > 

UM  •  delta 

•  PARAMETER  VARIATIONS  * 

A1  -  MAl-T  *  AlIC 
A2  ■  MAZ'T  ♦  A2IC 
B  •  MB*T  ♦  BIC 
K  ■  MK'T  ♦  KIC 

■  DISTURBANCE  MODEL 

W  -  WSINE*SINIRSW(LTVWF .MWFREQ«T.MWFREQ>*T)  ♦  ... 
WPULSE«PULSE(0. .UPULSP.WPULSW)  ♦  ... 

WNOISE*OU(RSW( LTVWNS. l/(MNOISE-T*.000l > . 1/NOISEW) , . . . 
WMEAN.WSTDEV)  ♦  ... 

WCONST 

•  PLANT  MODEL  ’ 

XID  ■  X2 

XI  ■  INTEGtXlD.Xl IC) 

X2D  -  -Al'Xl  -  A2-X2  ♦  U  <■  W 
X2  •  INTEG(X2D .X2IC  ) 

Y  -  B-Xl  *  K*X2 

•  REFERENCE  MODEL 


A-5 


XIMD  ■  X2M 


-  INTEG(XIMO.XIMIC) 

*.  -  IlI 

X2MD  -  -AN1»X1M  -  AN2*X2M  ♦  UM 
X2M  -  INTEG(X2MD,X2MIC) 

YM  -  BN*X1M  ♦  tCN*X2M 

•  ADAPTIVE  CONTROLLER 

EY  -  YM  -  Y 

CIO  -  TN*RSW(LERROR.EY. 1 . )*EY 
C20  ■  TN*RSW(LERROR.EY.l .  )«X1M 
C3D  ■  TN«RSW(LERROR,EY. 1 . )*X2M 
C40  -  TN«RSW(LERROR,EY, 1 . )«UM 

Ct  -  INTEGICIO.CIIC) 

C2  «  INTEG(C2D.C21C) 

C3  -  INTEG(C30.C3IC) 

C4  -  1NTEG(C4D.C41C) 

KE  -  TNB*RSW(LERROR.EY. 1 . >«EY  *  Cl 
KCl  -  TNB*RSU(LERR0R.EY.1. )»XIM  ♦  C2 
KC2  ■  TNB«RSy(LERROR.EY.l. )»X2M  ♦  C3 
KU  •  TNB«RSW(LERROR.EY.l. )«UM  ♦  C4 

U  -  RSW(LADAPT.RSW(LERRCH.XE«EY.j0r.  >.  . . 

*  RSWtLXMOD.KCl*XlM  ♦  XC2*X2M,« . > . . . 

♦  RSW(LUM.KU*UM.«. ) . UM ) 

•  TERMINATION 

ISE  ■  INTEG(EY«*2.«. ) 

TERMT  (dSE  .GE.  ISEMAX)  .OR.  (ABSIEY)  .GE.  EMAX)  ... 
•OR.  (T  .GE.  TSTOP)) 

END  S'OF  derivative* 

END  S'OF  dynamic* 

END 


A-6 


A. 4  PID  Simulation 


PROGRAM 

INITIAL 

'  ARRAY  QN(2) .00(2) 

LOGICAL  LDAC.LPULSE.LAOAC.LTVWF 
CONSTANT  AIIC-1.A2IC-2.AN1-1 .AN2-2 
CONSTANT  KIC-.S,BIC-1 .KN-.5.8N-1 
CONSTANT  MAI  -0.  .MA2-^.  .MB-JB.  .MK-flT. 

CONSTANT  KP-1..  RI-l.l.  KF-.55.  KO-BT. 

CONSTANT  XlIC-iB.,  X2IC-B. 

CONSTANT  XMlIC-flf.,  XM2IC-0. 

CONSTANT  CHIIC-B..  E2HIC-0. 

CONSTANT  ZAllC-Z..  ZA2IC-0. 

CONSTANT  21IC-0..  22IC-0. 

CONSTANT  ON-1 . .0. 

CONSTANT  OD-0. 1.1. 

CONSTANT  PULSEP«20. 

CONSTANT  TSTOP  -79 . 99 . ISEMAX- 1 . 3 . EMAX- . 1 
CONSTANT  WSINE-0. .WPULSE-0. ,WNOISE-0. .WCONST-0. 
CONSTANT  MWFREQ-0. 1 .WPULSP-40. .WPULSW-20. 

CONSTANT  WMEAN-0. .WSTOEV-l , .NOISEW-3. 

CONSTANT  L0AC-. FALSE. . L PULSE- .FALSE . 
constant  LA0AC-. false.  .  LTWF -.  FALSE  . 

PULSEW-PULSEP/2. 

U  -  0. 

CHl-0. 

END  S 

DYNAMIC 

DERIVATIVE 

CINTERVAL  CINT  -  0.05 
MAXTERVAL  MAXT  •  0.01 
NSTEPS  -  1 

■  INPUT  ■ 

DELTA  -  RSWC  LPULSE.PULSE(0. .PULSEP .PULSEW) , I . ) 

■  PARAMETER  VARIATIONS  * 

A1  -  MAl-T  ♦  AlIC 
A2  -  MA2*T  ♦  A2IC 
8  •  MB-T  ♦  BIC 

<  •  MK-T  ♦  KIC 

DISTURBANCE  MODEL 

W  -  WSlNE»SIN{RSW(LTWF.MWFREO*T.MWFREQ)*T)  ♦  ... 

»  WPULSE*PULSE(0. .WPULSP.UPULSW)  ♦  ... 

WNOISE«OU( ( 1/NOISEW) .WMEAN.WSTDEV)  ♦  ... 

WCONST 

■  CONTROL 

E  -  DELTA  -  Y 

U  -  KF-(kP-E  ♦  KI-INTEG(E.0. )  ♦  KD*TRAN ( 1 . 1 . ON . OD . E ) > 
•  PLANT  MODEL 

XID  -  X2 

XI  -  INTEG(XID.XIIC) 


A-/ 


X20  -  -Al*Xl  -  A2«X2  ♦  U  ♦  W 
X2  -  INTefiCX20,X2IC) 

Y  -  K*{ (B/K)«X1  ♦  X2) 

•  IDEAL  MODEL  * 

XMl^-  INTEGIXMIO.XMI IC> 

XM2D  -  -AN1«XMI  -  AN2*XM2  ♦  DELTA 
XM2  -  1NTEG(XM2D.XM21C) 

YM  -  <N*{ (BH/KN)«XM1  ♦  XM2) 

•  INTEGRAL  SOUAREO  ERROR  ‘ 

ERROR  -  YM  -  Y 

ISE  -  INTEG(ERROR**2.jB.  > 

•  termination 

TERMT  tdSE  .GE.  ISEMAX)  .OR,  (ABSIERROR)  .SE. 
•OR.  a  .GE.  TSTOP)) 


END  $*OF  DERIVATIVE* 
END  S’OF  DYNAMIC* 

END 


EMAX) 
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A. 5  STR  Simulation 


PROGRAM 

INITIAL 

ARRAY  QNUY(A).  QOUYIA).  QNYRtFI*).  QDYREFIA) 

LOGICAL  LPULSE.LSINE.LNOISE.LSTEP 
LOGICAL  DMOOEL .LCLOOP .TRUPAK.LGAIN 
ClNTERVAL  CINT  -  .5 
MAXTERVAL  MAXT  -  .Bl 
NSTEPS  NSTP  -  1 

constant  lambda  -  1  . 

CONSTANT  COVIC-lflfB. 

CONSTANT  TSTOP-24.99 

CONSTANT  TSAMP-.01 

CONSTANT  THTlIC-0..  THT2IC-#. 

CONSTANT  THT31C-.5.  THT4IC-1. 

CONSTANT  PULSEP-20. .PULSEW-10, ,FS1NE*.25 
CONSTANT  KPULSE-l .  .RSINE-l  . 

CONSTANT  RMSN-.5 

CONSTANT  A1 IC-1 . A2 I C -2 . AN  1  - 1 , AN2-2 
CONSTANT  KIC-.S.BIC-l.KN-.S.BN-l 
CONSTANT  MAl-0. .MA2-0. ,MB-0. 

CONSTANT  XllC-0..  X2IC-0. 

CONSTANT  XlMIC-0.,  X2MIC-0, 

CONSTANT  YnC-0..  Y2IC-0. 

CONSTANT  UlIC-0..  U2IC-0. 

CONSTANT  XlHIC-0.,  X2HIC-0. 

CONSTANT  21H1C-0..  Z2HIC-0..  23HIC-0. 

CONSTANT  W«0. 

CONSTANT  LPULSE«.TRUE. .  LS I NE« . FALSE . 

CONSTANT  LNOISE-. FALSE. .  LSTEP- . FALSE . 

CONSTANT  DMODEL-. FALSE.  .  LCLOOP- . FALSE TRUPAR- . TRUE . 
CONSTANT  LGAI N- . FALSE . 

PULSEW-PULSEP/2 . 

YL 1  -  0. 

Yl2  -  0. 

L'-l  -  0. 

UL2  -  0. 

UN  -  0. 

UF  .0, 

UCL  -  0. 

Uv  •  0. 

THTIHL  -  THTIIC 
THT2HL  ■  THT2IC 
THT3HL  •  THT3IC 
THT4HL  -  THT4IC 
PUL  ■  COVIC 
P22L  •  COVIC 
P33L  •  COVIC 
P44L  «  COVIC 
P12L  -  0. 

P13L  ■  0. 

P 1 4L  •  0. 

P23L  •  0. 

P24L  •  0. 

P34L  •  0. 

GSS-1  . 

KH-KN 
8H-BN 
A  1  H-AN  ; 

A2h-AN2 

END  $ 

DYNAMIC 
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DERIVATIVE 


•  PARAMETER  VARIATIONS  * 

Al  -  MAl«T  ♦  AlIC 
A2  -  MA2*T  ♦  A2IC 
B  ■  MB'T  ♦  BIC 
K  •  MK*T  ♦  KIC 


•  INPUT  • 


VREF  -  RSW(LPULSE.KPOLSE«(PULSE(0. .PULSEP.PULSEW)  > .  I  .  > 
RSW(LSINE.<SINE“SIN{2  =  *3.U15*FSINE“T),flr.)  *  ... 

RSW!LSTEP,1..0.) 

PROCEDURAL 

QNVREFC  n«l  . 

QNYREF(2)-A2H 
QNYREF ( 3 )»A1H 
QDYREFI 1  )-l  . 

QOYREP  (  2  )  -J?. 

QDYREF(3)«iB. 

END 

U  -  TRAN(2, 2. QNYREF. QOYREF.YREF)  -  RSWdCUOOP  .UY.flT.  ) 

•  PLANT  MODEL  * 

XID  -  X2 

XI  «  INTEG(XID.XIIC) 

X2D  ■  -Al'Xl  -  A2*X2  ♦  U  ♦  W 
X2  -  INTEG(X2D.X2IC) 

Y  -  B*Xl  ♦  K*X2 

•  filter  Y  for  use  in  ESTIMATOR  * 

YIO  -  Y2 

Yl  -  INTEG( YID . Yl IC ) 

Y2D  •  -AN2*Y1  -  AN1-Y2  ♦  Y 
Y2  -  INTEG{Y20,Y2IC) 

UID  -  U2 

U1  -  INTEGIUIO.UIIC) 

U2D  •  -AN2*U1  -  AN1-U2  *  U 
U2  -  INTEG(U20,U2IC) 

■  PLANT  MODEL  ’ 

XIMD  -  X2M 

XIM  -  INTEG(XIMO.XIMIC) 

X2M0  -  -AN1«X1M  -  AN2*X2M  *  YREF 
X2M  -  INTEG(X2MD,X2MIC> 

YM  -  BN*X1M  ♦  XN«X2M 

■  CONTROL  LAW  “ 

PROCEDURAL 

QNUY( 1 >-2. 

QNUY(2)-1 .♦2«A2H 


:  if  QNUY( 3>-A2*2»AlH 
■  —  -  ■  QNUY(4>-AIH 
ODUY{ 1 )-KH 
QDUY( 2 )-BH 
QDUY(3)-^. 

QDUY( 4 )-0. 

END 

UY  ■  TRAN{3.3.QNUY,Q0UY.Y) 

END  S*OF  DERIVATIVE* 

DISCRETE 

INTERVAL  TS-.01 
PROCEDURAL 

•  PARAMETER  ESTIMATION  USING  RECURSIVE  LEAST  SQUARES  ... 
WITH  EXPONENTIAL  FORGETTING  * 

E  -  Y  -  ( -THT1HL*Y2-THT2HL«Y1+THT3HL»U2*THT4HL*U1 ) 

W1  -  -PI 1L»Y2-P12L»Y1*P13L*U2*P14L*U1 

W2  »  -P 12L«Y2-P22L«Y1*P23L*U2+P24L*UI 

W3  -  -P13L*Y2-P23L*Y1*P33L*U2*P34L*U1 

W4  -  -P14L''Y2~P24L»Y1+P34L*U2*P44L''U1 

DEN  •  -W1<'Y2-W2*Y1*W3*U2*W4»U1*LAM80A 

THTIHN  -  THTIHL  ♦  W1*£/0EN 
•.■HT2HN  -  THT2HL  *  W2»E/DEN 
THT3HN  -  THT3HL  ♦  W3*E/DEN 
THT4HN  •  THT4HL  ♦  W4-E/DEN 

PUN  -  (PllL  -  (Wl**2)/DEN)/LAMB0A 
P22N  •  (P22L  -  (W2*«2 )/OEN)/LAMgOA 
P33N  -  (P33L  -  {W3*«2 )/DEN )/LAMB0A 
P44N  -  (P44L  -  (W4«»2 )/DEN)/LAMBDA 
P12N  •  <P12L  -  (W1«W2)/DEN)/LAMBDA 
P13N  •  (P13L  -  (W1«W3)/DEN)/LAMB0A 
PUN  -  (PUL  -  (W1«W4)/DEN)/LAMB0A 
P23N  -  (P23L  -  (W3»W2 ) /DEN >/LAMBDA 
P24N  -  (P24L  -  (W4»U2 )/DCN >/LAMBDA 
P34N  ■  (P34L  -  (W3*U4)/DEN)/LAMRDA 

A2H  -  ANl  ♦  THTIHL 
AlH  -  AN2  ♦  THT2HL 
XH  ■  THT3HL 
BH  -  THT4HL 

THTIHL  ■  THTIHN 
THT2HL  -  THT2HN 
THT3HL  ■  THT3HN 
THT4HL  ■  THT4HN 
PUL  ■  PUN 
P22L  ■  P22N 
P33L  -  P33N 
P44L  •  P44N 
P12L  ■  P12N 
P13L  -  P13N 
PUL  -  PUN 
P23L  -  P23N 
P24L  -  P24N 
P34L  ■  P34N 

END  S*OF  PROCEDURAL" 


END  »'0F  DISCRETE* 

•  TERMINATION 

TERMT  (T  .GE.  TSTOP ) 
END  S*OF  DYNAMIC* 

END 
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